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Quenching of Leading Jets and Particles: the p^ Dependent 
Landau-Pomeranchuk-Migdal effect from Nonhnear k± Factorization 
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We report the first derivation of radiative nuclear stopping (tiie Landau-Pomerancliuk-Migdal 
effect) for leading jets at fixed values of the transverse momentum p in the beam fragmentation 
region of hadron- nucleus collisions from RHIC (Relativistic Heavy Ion Collider) to LHC (Large 
Hadron Collider). The major novelty of this work is a derivation of the missing virtual radiative 
pQCD correction to these processes - the real-emission radiative corrections are already available in 
the literature. We manifestly implement the unitarity relation, which in the simplest form requires 
that upon summing over the virtual and real-emission corrections the total number of scattered 
quarks must exactly equal unity. For the free-nucleon target, the leading jet spectrum is shown to 
satisfy the familiar linear Balitsky-Fadin-Kuraev-Lipatov leading log-i (LL-) evolution. For nuclear 
targets, the nonlinear fcx-factorization for the LL-i evolution of the leading jet spectrum is shown 
to exactly match the equally nonlinear LL- evolution of the collective nuclear glue - there emerges 
a unique linear fc_L-factorization relation between the two nonlinear evolving nuclear observables. 
We argue that within the standard dilute uncorrelated nucleonic gas treatment of heavy nuclei, in 
the finite energy range from RHIC to LHC, the leading jet spectrum can be evolved in the LLi 
Balitsky-Kovchegov approximation. We comment on the extension of these results to, and their 
possible reggeon field theory interpretation for, mid-rapidity jets. 



PACS numbers: 13.87.-a, 11.80La,12.38.Bx, 13.85.-t 



I. INTRODUCTION 

In this communication we address two related 
issues in the perturbative Quantum Chromo Dy- 
namics (pQCD) description of the single-jet in- 
clusive spectra in hadron-nucleon and hadron- 
nucleus collisions. Our principal interest is in the 
nonlinear fc_L factorization for the pQCD radia- 
tive quenching (stopping) of leading jets with fixed 
transverse momentum p produced at large values 
of the Feynman variable xp, i-e., in the beam frag- 
mentation region of hadron-nucleon and hadron- 
nucleus collisions. As it has become customary 
for the nuclear dependence of radiative effects, we 
refer to nuclear stopping of leading jets at fixed p 
as the p-dependent Landau-Pomeranchuk-Migdal 
(LPM) effect. There is an extensive literature on 
the LPM effect for mid-rapidity jets, where the 
prime concern is in the radiation induced loss of 
the transverse, with respect to the hadron-nucleus 
collision axis, momentum of the jet when one in- 
tegrates over the transverse momentum of the ra- 
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diation with respect to the jet axis (0], for re- 
views and further references see Q, |^). In the 
mid-rapidity kinematics, the transverse momen- 
tum of radiation with respect to the parent high- 
p±^ parton only amounts to a small shift of the 
(pseudo)rapidity of the high-pj^ parton which is 
hardly observable due to the boost invariance in 
the mid-rapidity region. In contrast to that, lead- 
ing jets with xp ~ 1 are generated by hard scat- 
tering of valence quarks. The density of valence 
quarks drops rapidly at large values of the Bjorken 
variable in the beam nucleon, xn — *• 1, which en- 
tails that the xp dependence of the leading jet 
cross section is extremely sensitive to the radiative 
energy loss. We develop a nonlinear k^ factoriza- 
tion description of leading jet production with full 
allowance for the real and virtual QCD radiative 
corrections. In the second part we focus on the 
leading log^ (LL-) properties of the leading jet 
spectrum. We argue that this spectrum offers a 
long sought linear mapping of the unintegrated 
collective nuclear glue and comment on the ex- 
tension of this finding to mid-rapidity jets. 

Within pQCD, there are several subtleties in 
the application of factorization theorems to the 
leading quark production in qN, qA collisions. 
In the more familiar coUinear approximation de- 



scribed in all textbooks |j|, the relevant pQCD 
subprocess is hard scattering of valence quarks 
off soft gluons with x ^ 1 in the target, which 
can also be viewed as a hard breakup (excita- 
tion) q — > q'g. All partons enter the hard pQCD 
scattering with vanishing transverse momentum 
and the large p^ of the observed jet is gener- 
ated in the real-emission hard breakup. The 
inadequacy of such a coUinear approximation is 
well known, a recent concise discussion of the 
necessity of fully unintegrated parton densities 
with explicit allowance for transverse momenta is 
found in p. In the high-energy hmit, the lat- 
ter is furnished by the so-called fc^-factorization 
which goes back to the classic 1975-1978 works 
on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) 



equation for the small- a; evolution of the unin- 
tegrated glue (la, Q; see also the recent review 
on fc^-factorization and more references in [al3)- 
In contrast to the coUinear factorization, within 
the fc^-factorization, the hard breakup q — > q'g 
is but the real-emission pQCD radiative correc- 
tion (REC) to the lower order (Born) subprocess, 
which is the quasielastic scattering qN — > qN* , 
qA -^ qA* , where the target debris N* , A* is 
left in a color excited state. Evidently, the real- 
emission must be complemented by the match- 
ing virtual pQCD radiative corrections (VRC) to 
the lower order quasielastic scattering - the real- 
emission and virtual radiative corrections are two 
integral parts of the small-x evolution of quasielas- 
tic scattering, see Fig. ^ 
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FIG. 1: The real emission (REC) and virtual radiation (VRC) contributions to the Born spectrum of leading 
quark jets in inclusive production off nuclei. The (pseudo) rapidity of radiation correction gluons runs between 
the quark jet rapidity and the nuclear boundary condition rapidity tia = log — , where xa is defined by Eq. 
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Hard scattering off free nucleons is described 
by single-gluon exchange in the t-cannel. In colli- 
sions with nuclei, the enhancement of multigiuon 
exchanges leads to a novel co ncep t of nonlinear 
fcj^ factorization [Ulilllllillllli^. The Born 
and VRC amplitudes describe final states in which 
the large transverse momentum of the leading jet 
is compensated by hadrons at rapidities r] < rjA, 
while in REC such a compensation is provided 
also by gluon jets with rapidities "Ha "^"Hg '^ Vq' ■ 
A theoretical treatment of the p dependent LPM 
effect for leading jets is necessary for the interpre- 



tation of the forthcoming experimental data on 
leading particles from the RHIC to LHC ([ig and 
references therein). 

The interplay of the virtual and real radiative 
corrections to single jet spectra is controlled by 
the unitarity condition. In its simplest form uni- 
tarity demands that the multiplicity of final state 
leading quarks from the real radiation, and from 
the VRC-corrected lower order (Born) quasielas- 
tic scattering add up exactly to unity. In other 
words, the small- a; evolutions of the total inclu- 
sive cross section for the leading quark production 



and of the total cross section of the quark-target 
interaction must preserve their exact equahty. Ev- 
idently, the accurate evaluation of the energy loss 
requires the REC for a finite fraction, Zg, of the in- 
cident quark's momentum carried by the radiated 
gluon. Correspondingly, the VRC must be calcu- 
lated to a matching accuracy. While the required 
real-emission spectra do exist in the literature for 
both the free-nucleon and nuclear targets (we fol- 
low the approach of Ref. iTnl . fo r related works 
and further references see [H, IS 123 ) j the match- 
ing derivation of virtual radiative corrections was 
not yet available - it is a major novelty of this 
paper. Our derivation of the combined VRC and 
REC fully respects the s-channel unitarity con- 
straints and paves the way to a first consistent 
description of the LPM effect for leading jets for 
fixed p of the jet. 

Much insight into the interplay of the virtual 
and real-emission radiative corrections for nuclear 
targets comes from the study of LLi evolution of 
the LPM effect-deconvoluted inclusive spectra of 
leading particles, i.e., neglecting the radiative en- 
ergy loss. At moderately high energies, the unitar- 
ity constraints in hard qN collisions are still unim- 
portant and they can be described by single-gluon 
exchange in the i-channel between the projectile 
and target. In this case a comparison of the LL- 
evolution of the LPM-effect deconvoluted leading 
parton spectrum and of the BFKL evolution of the 
unintegrated glue of the nucleon leads to a simple 
fc^-factorization relationship 



^^^ = - fix p) 
d^p 2^^ '^' 



(1) 



(for the definition of the unintegrated glue f{x,p) 
in the free nucleon see below). I.e., the leading 
quark spectrum satisfies the same BFKL small- a; 
evolution equation as the unintegrated glue [3, Q , 
and this spectrum emerges as a direct probe of 
the unintegrated glue in the target nucleon. It 
is remarkable that our derivation of this prop- 
erty makes a direct use the s-channel unitar- 
ity for color dipole scattering. In the realm of 
the BFKL approach, this property can be re- 
garded as a straightforward consequence of the 
correspondence between LL^ approximation and 
the dominance of multiregge production processes 

IEIEl2i|. 

In contrast to the single t-channel gluon ex- 
change with the free-nucleon target, multiple 
gluon exchanges with the target nucleus are en- 
hanced by a large thickness of the nucleus. Here 
the very possibility of describing a nucleus by a 
unique collective unintegrated glue, and the ex- 



istence of factorization theorems, become prob- 
lematic. Various aspects of this issue have been 
addressed to within the so-called Color Glass 
Condensate approach, in which the starting idea 
was to describe the nucleus by the collective 
Weizsacker- Williams gluon density ( |22| and refer- 
ences therein). More recently, much progress has 
been made within the nonlinear fcj^-factorization 
reformulation of the multiple scattering theory for 
color dipoles. Within the color dipole approach, 
the extensive previous studies of the dijet spectra 
(Up, in El 111 II1II3, see also the recent Ref. 
|23j) and single-jet spectra [i3 revealed a striking 
breaking of the conventional linear fcj^ factoriza- 
tion for hard processes in a nuclear environment. 
Specifically, all the dijet observables can be repre- 
sented as nonlinear quadratures of the collective 
nuclear glue defined in terms of a reference nuclear 

§cess - production of coherent diffractive d i.jets 
m. On the other hand, in Refs. [ll[l3 
we noticed the phenomenon of Abelianization of 
the single particle spectra - to LL^ the nonlin- 
ear fc_L-factorization simplifies down to the linear 
fcj^-factorization subject to the judicious choice of 
the component of the color space density matrix 
for the collective nuclear glue. Here we demon- 
strate that, at least to the first order of LL^, a 
similar Abelianization takes place for the leading 
jet spectrum from inelastic processes. 

When speaking of fc^-factorization for nuclear 
collisions, one needs to define the collective nu- 
clear glue in terms of the properly chosen physics 
observable. Here we start with the BFKL defi- 
nition of the unintegrated glue for free nucleons 
la, • As noticed in J2g , this unintegrated glue 
can be directly measured in diffractive produc- 
tion of dijets (for other applications of the un- 
integrated glue see the reviews in 0, I3, H, 123). 
Generalizing this observation to nuclear targets 
\WL |2J|, we define the collective nuclear glue in 
terms of the color-dipole nuclear S-matrix which 
is accessible experimentally via diffractive dijets. 
The LLi evolution of the so-defined glue is a 
nonlinear one. The REC to the inclusive sin- 
gle leading parton spectrum derived in [l7| is a 
nonlinear functional of the collective nuclear glue. 
Equally highly nonlinear is the VRC correction to 
quasielastic parton-nucleus scattering derived in 
this paper. Our nontrivial finding is that the lead- 
ing parton spectrum and the collective glue share 
identical first iteration of the nonlinear LL- evo- 

X 

lution, which is the second novelty of our paper. 
One is tempted to conjecture that the spectrum 
of leading jets is linear fc_L-factorizable in terms 
of the collective nuclear glue to all order of the 



nonlinear LL-!- evolution, i.e., it is a long sought 
unique linear probe of the collective nuclear glue. 
It is by now understood that the closed form of 
all-order LL^ evolution for nuclei does not exist 

X 

(a detailed review is found in ^J ) . Still we argue 
that in view of a high accuracy of the dilute un- 
correlated nucleonic gas picture of heavy nuclei, 
the small- a; evolution of the leading jet spectrum 
from RHIC to LHC can to a controlled accuracy 
be described by several iterations of the Balitsky- 
Kovchegov (BK) approximation 29, 30] . 

We find it very instructive to compare nonlin- 
ear fc_L-factorization properties of the LL- evolu- 
tion of the inclusive spectrum of leading quarks 
belonging to the color-triplet (fundamental) and 
color-octet (adjoint) representations of SU{Nc). 
The latter case is of interest from the viewpoint 
of supersymmetric QCD, one can also think of 
leading gluons in interactions of glueballs. A pat- 
tern of the nonlinear fc_L -factorization for the spec- 
trum changes strikingly from the triplet to octet 
quarks, equally different is the nonlinear BFKL 
evolution of the components of the color den- 
sity matrix for the collective glue for the triplet- 
antitriplet and octet-octet color dipoles. For in- 
stance, for strongly absorbing nucleus and triplet- 
antitriplet dipoles, the nonlinear BFKL evolution 
looks as a fusion of two nuclear pomerons, i.e., 
it is dominated by triple-pomeron interactions. 
In striking contrast to that, for the octet-octet 
dipoles the dominant term looks as a fusion of 
three pomerons from one component of the color 
density matrix for the collective nuclear glue to 
the fourth pomeron described by a different col- 
lective glue. In both cases, however, the evolution 
of the two quantities - the leading parton spec- 
trum and the collective glue - turns out to be 
identical, which lends more support to our con- 
jecture on leading jets as a linear probe of the 
collective nuclear glue. When reinterpreted in the 
reggeon field theory language, our results corre- 
spond to a full resummation of multiple pomeron 
exchanges enhanced by a large thickness of the 
target nucleus. Such a resummation is interesting 
on its own, an extension of these results beyond 
the dilute projectile and the Balitsky-Kovchegov 
approximations, and the challenging inclusion of 
the pomeron loop effects (for the review see [23 |') 
go beyond the scope of the present communica- 
tion. 

The presentation of the main material is orga- 
nized as follows. In Sec. II we introduce the basic 
color dipole formalism, a definition of the collec- 
tive nuclear unintegrated glue and derive the dif- 
ferential cross section of quasielastic scattering of 



quarks off free nucleons and nuclei to the Born 
approximation. The color-dipole master formulas 
for the real and virtual radiative corrections to the 
inclusive spectrum of leading quarks are presented 
in Sec. III. Our derivation of the virtual correction 
in Sec. III.B makes a manifest use of the unitarity 
relation. We then proceed to the (nonlinear) k±- 
factorization for the spectrum of leading quarks 
off the free nucleon (Sec. IV. A) and color-triplet 
quark-nucleus (Sec. IV. B) collisions as well as for 
the spectrum of leading color-octet quarks (and 
gluons) in collisions with nuclei (Sec. IV. C). These 
results arc sufficient for the quantitative treatment 
of the LPM effect at fixed transverse momentum 
of leading jets at RHIC energies. In Sec. IV. D 
we revisit the issue of unitarity and demonstrate 
the consistency of our radiative corrections with 
the leading quark number sum rule. Inspired by 
the observations in Sec. II and Sec. IV, in Sec. V 
we focus on the LLj evolution properties of inclu- 
sive single particle spectra, which are necessary for 
the extension of our technique to LHC energies. 
Wc demonstrate, that to the LLj approximation, 
the nonlinear-evolving leading parton spectrum in 
parton-nucleus collisions proves to be linear k±- 
factorizable in terms of the nonlinear-evolved col- 
lective nuclear glue. In Sec. V.D we comment 
on the unified description of the LL- evolution of 
the triplet-antitriplet and octet-octet S-matrices. 
The subject of section V.E is a partial vindica- 
tion of the BK approximation - while it is not a 
closed-form equation and is not really applicable 
to the free- nucleon target, for heavy nuclei sev- 
eral iterations of the BK approximation will be 
a useful working approximation in the range of 
energies from RHIC to LHC. This gives a practi- 
cal procedure for the quantitative evaluations of 
the LPM effect for leading jets at LHC. In section 
V.F we comment on how our treatment of the in- 
terplay between VRC and RFC extends to the 
production of mid-rapidity jets, and on the possi- 
ble interpretation of our results from the effective 
reggeon field theory perspective. In the Conclu- 
sions we summarize our principal findings. In the 
Appendix we comment on an interesting observa- 
tion that the nonlinear fc_L -factorization compo- 
nent of the LLi nonlinear evolution gives a pure 
higher twist contribution to the collective nuclear 
glue 113. 
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FIG. 2: (a,b) The Feynman diagrams for Upsilon 
breakup by single-gluon exchange with the target nu- 
cleon and (c) an archetype diagram of breakup of Up- 
silon by a multiple-gluon exchange with the target nu- 
cleus. 



II. LEADING PARTICLE SPECTRA TO 
THE BORN APPROXIMATION 

A. Leading particles from the breakup of 
heavy quarkonium 

To set up the framework, we start with a toy 
problem: the breakup of the color-singlet heavy 
quarkonium (let it be T) in a high energy col- 
lision with the free-nucleon target TA^ -^ bbX. 
Heavy quarkonium can well be approximated by 
its quark-antiquark Fock state with the size Rr ^ 
1/mf,. The lowest order pQCD subprocess is a 
breakup by Upsilon-gluon fusion, Tg't — *■ bb, where 
gt is the gluon exchanged in the i-channel. (Fig. 
|21l. The corresponding inclusive differential cross 
section, obtained after summing over all excita- 
tions of the target nucleon, is well known; it is 
found, for instance, in ^^ (Eqs. (4) and (5)): 



dzd^pcFA 2TrNc A^ 

x|*t(^,p)-*t(2,p- A)|2. 



(2) 



Here asip^) is the pQCD coupling, Nc is the num- 
ber of colors, ^t(-z,p) is the lightcone wave func- 
tion of the quarkonium, p is the transverse mo- 
mentum of the quark in the quarkonium, z is the 
fraction of the lightcone momentum of the quarko- 
nium carried by the quark, A is the total trans- 
verse momentum of the bb pair. 



Ti^x, K ) 



2, _ dG{x,K^^ 



d log K^ 



(3) 



is the so-called unintegrated gluon density in the 
target nucleon as introduced in the BFKL ap- 
proach ja, Q- The Bjorken variable x of the ex- 



changed gluon equals 



^ll 



M| 



W^lp 



(4) 



where PFxp is the cms energy in the Tp system and 
Mf^l is the transverse mass of the excited quark- 
antiquark pair. 



Mil 



iP 



(5) 



The transverse momentum A of the bb pair comes 
entirely from the exchanged gluon and the differ- 
ential cross section Q emerges as a natural probe 
of the unintegrated glue J^{x, A^). We recall that 
for the free nucleon target J-{x,k^) is a solution 
of the LLi BFKL evolution equation 0,13' ^ ^^~ 
tailed discussion of the connection between the in- 
tegrated glue G'(x, K^) of Eq. Q and the gluon 
densities of the leading order Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi evolution equations J32| is 
found in Refs. [liilllllllli and need not be 
repeated here. 

Now consider transverse momenta p much 
larger than the typical momentum of the Fermi 
motion in the quarkonium, pf, ~ l/i?x- The wave 
function ^(z,p) vanishes steeply beyond pf,, and 
the single-quark spectrum takes the convolution 
form 



da{T 



bb) f ,2 .^s(p^) T{x,K?) 



^ / d^K- 



dzd^p J " 2TrNc 

x\-^riz,p)-'^T{z,p-K)\^ 



d^K 



as[p^) T{x,K?) 



2ttNc /«4 

/ d^Kd^Pf,S{p — p^ — K 
dqb/r{z,Pb) da{b 



*t(2,p-k)|^ 



b') 



Here 



d^Pb 



dqb/riz,Pb) 
d^Pb 



d'^t 



= l*T(z,Pfc)| 



(6) 



(7) 



is the momentum distribution of quarks in the 
quarkonium and 



d(T{b-^b') as{p^) T{x,K^ 



d^t 



2llNr 



(8) 



must be interpreted as a differential cross sec- 
tion of the valence quark quasielastic scattering 
bN -^ b'X. Furthermore, the longitudinal mo- 
mentum distribution of leading quarks is given by 



precisely the transverse momentum-dependent z- 
distribution in the quarkonium - z is the Feyn- 
man variable of the observed quark. As such, the 
convolution (|SJ) can be regarded as the archetype 
fcj^-factorization for the production of hard lead- 
ing quark jets in high-energy hadron-nucleon col- 
lisions. The adjective 'hard' refers to the jet mo- 
menta harder than the intrinsic momentum of 
quarks in the quarkonium: p^ ^ {pI)- Now we 
proceed to the fcj^-factorization properties of the 
quasielastic scattering of partons - quarks and 
gluons - off free nucleons and heavy nuclei. 



B. The free-nucleon target: the dipole cross 

section, unintegrated glue and quasielastic aN 

scattering 

We start with the dipole cross section for the 
aa color dipole. It is described by two-gluon ex- 
change in the i-channel (Fig. |3). The S-matrices 
of the aN and aN interaction at impact parameter 
ba{ba) equal, respectively [W^ . 



Siib,) = l-zT:V^xib-a)-lT':T':x'ib,), 



(9) 



were T"Vax{b) is the eikonal operator for the aN 
single-gluon exchange interaction. In the realm 
of the Standard Model QCD, quarks belong to 
the color-triplet (fundamental) representation of 
SU{Nc), but much hght on the small- a; evolution 
properties of the nuclear spectra is shed by consid- 
ering quarks also in the octet (adjoint) represen- 
tation of SU{Nc). One can also think of leading 
gluons in interactions of glueballs. The expansion 
lO to the considered pQCD approximation satis- 
fies the unitarity condition 



Saiba)SUba) = t. 



(10) 



The spin of partons is unimportant in view of 
the s-channel helicity conservation in high energy 
QCD. The important observation is that S'^ equals 
the S-matrix for interaction of the antiparticle a 
[Sft ISTj . The vertex operator Va for excitation of 
the nucleon g^N -^ N^ into a color octet state N* 
is so normalized that after the application of clo- 
sure over the final state excitations TV* the vertex 
g°'gl^NN equals {N\V^Vf3\N) = S^p. The second 
order terms in do already use this normaliza- 
tion. 
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FIG. 3: The four Feynman diagrams for the quark- 
antiquark dipole-nucleon interaction by the two-gluon 
pomeron exchange in the i-channel. 



The S-matrix of the (aa)-nucleon interaction 
equals 



Saci{ba,ba) 



_ {N\TT[Sa{ba)Si{b-a)]\N) 
(Af|Trl|iV) 

l~\Ca[x{ba)~X{b-a)?. (H) 



where Ca is the quadratic Casimir operator. The 
corresponding profile function is Taa{ba,ba) — 
1 — Saa{ba^ ba), and the dipole cross section for in- 
teraction of the color-singlet aa dipole r = ba~ba 
with the free nucleon is obtained upon the integra- 
tion over the overall impact parameter 



(Jaa[x,r 



2 j SbTaa{b,b-r) 

Cajd?b[x{b)-x{b-r)f 



= 7^cr(a;,r), 



(12) 



where Cp = {N^ — i)/2Nc is the quark Casimir 
and a{x,r) is the dipole cross section for the qq- 
dipole. It is related to the gluon density in the 
target by the fc^-factorization formula |3J, [Sg 

a{x,r) = cpKf{x,K.)[l — cxp{iKr)] , (13) 

where 

/(:r,K) = ^.1 T{x,k'). (14) 

Nc K^ 

The leading Log^ evolution of the dipole cross 
section is gove rned by the color-dipole BFKL evo- 
lution 13^ Ispl , the same evolution for the uninte- 
grated gluon density is governed by the familiar 
momentum-s] 
large dipoles 



momentum-space BFKL equation 



by the tan 
B,U For 



<ya-aA^) = CaJ d^b[x\b)+X^{b^r)] 



very 



- 2Ca / d?b 



(A^|Tr{a-Sa(b)}|iV} 



(A^|Trl|A^) 



C„ 



2aa{x) = -^oo{x), 



(15) 



where cro{x) is the cross section for large qq 



dipoles. 

Now we turn to the master formula for the inclu- 
sive differential cross section of quasielastic scat- 
tering aN -^ a'N*. Only the final states N* ^ N 
are excited by the f-channel gluon exchange and 



dar 



d?p 



{2^Y 



(fb(fb'eyip[-~ip{b 






(7V|Tr{[a - S,(b)]|7V*)(7Vl[l ~ Si (6')]}|A^) 



(A^lTrUlTV) 



(16) 



r 



After the application of closure 



^ \N*){N* 



\N){N\^l, 



(17) 



and upon simple algebra one finds 



daQei{x) 



—-^ I d2cexp[-ipc] 



d2p 2(27r)2 

X [(Jaafl{x) -(Taa{x,c)] 
Ca 



2Cf 



f{x,p), 



(18) 



where p is the transverse momentum of the scat- 
tered quark, cf. with Eq. ©. To this Born 
approximation, quasielastic scattering off the nu- 
cleon exhausts the total cross section, the pure 
elastic cross section is of higher order in pQCD 
perturbation theory - it starts with four gluons in 
the t-channel. Quasielastic qN scattering emerges 
as a direct probe of the unintegrated gluon struc- 
ture function of the target nucleon ~ we shall show 
this to hold beyond the Born approximation. The 
above explicit calculation of the emerging single- 
nucleon matrix elements is straightforward and 
hereafter, unless that might cause a confusion, it 
will be suppressed. 



C. Quasielastic qA scattering and collective 
nuclear glue 

Here we need to specify the boundary condition 
of the small-x evolution. In the interaction with 
heavy nuclei the coherent nuclear effects evolve, 
and the concept of the collective unintegrated nu- 
clear glue becomes applicable, for partons with 



X <XA 



1 



2i?yiTOjv 



(19) 



In the Breit frame, this coherence condition cor- 
responds to the spatial overlap and fusion of par- 
tons belonging to different nucleons at the same 
impact parameter in the Lorentz-contracted ultra- 
relativistic nucleus |40j . At a; > x^ coherent nu- 
clear effects are absent and the impulse approx- 
imation holds. The small- a:: evolution down to 
XA, the span of which could be substantial for 
the parametrically small xa for very heavy nuclei, 
must be performed at the free nucleon level. The 
evaluation of coherent nuclear effects at x < xa 
must start with the boundary conditions set by 
the free- nucleon quantities evaluated at xq = xa- 
We refer to this boundary condition as the nuclear 
Born approximation. 

Interaction of the generic n-parton system with 
the target nucleon is described by the coupled- 
channel cross section operator E^"^. In the regime 
of coherent interaction with the heavy nucleus, 
X < Xa, the corresponding S-matrix operator is 



given by the Glauber-Gribov formula 



operato 



where 



T{b) 



dr^UAib^rz) 



(20) 



(21) 



is the optical thickness of the nucleus. The nu- 
clear density nA{b, rz) is normalized according to 

d^fuAib, Tz) = I d^bT{b) = A, where A is the 



nuclear mass number. The principal assumption 
behind H2U|I is that the nucleus is a dilute uncor- 
related gas of color-singlet nucleons. Also, color 
dipoles are much smaller than the size of the nu- 
cleus. All the nuclear cross sections are defined 
per unit area in the impact parameter space. 



For the transformation from the color-dipole to 
the momentum representation, we define the col- 
lective nuclear unintegrated gluon density per unit 
area in the impact parameter plane, (/i(6,xo, /«), 
in terms of the Born-approximation nuclear pro- 
file function for the triplet-antitriplet qq dipole 

EaHISIli: 



r[6, cr(xo,r)] = 1 -S[5;cr(a;o,r) 



= 1 — cxp 



--aixo,r)T{b) 



d^K(j){b, xq,k) 1 — exp{iKr 



(22) 

The expansion for (j){b,xo,K) in terms of the 
collective glue for overlapping nucleons in the 
Lorentz-contracted ultrarelativistic nucleus, and 
its nuclear shadowing and antishadowing proper- 
ties, are found in [i3,ll3- The utility of 0(6, xq, k) 
stems from the observation that the driving term 
of small- a: nuclear structure functions, the am- 
plitude of coherent diffractive production of di- 
jets off nuclei and the single-quark spectrum from 
the 7* — >• qq excitation off a nucleus all take the 
familiar linear fc_L-factorization form in terms of 



(j)(b,xo, k). Still more convenient, and more uni- 
versal, would be a definition in terms of the color- 
dipole S matrix, 

$(6,xo,/«) = / S[b, g(a;o,'r)]exp(-tAt-r) 

= S[h,aQ{xo)]5^'^\K,) + (l){h,xo,K), 

(23) 



which by its definition satisfies the sum rule 



/ (fK(^{b,xo,n) = 1. 



(24) 



Now we proceed to the differential cross section 
of quasielastic scattering a A —f a A* . Heavy nu- 
clei are strongly absorbing targets and the optical 
nuclear thickness T(b) is a new large parameter 
in the problem. Because of enhanced multiple t- 
channel gluon exchanges the pure elastic scatter- 
ing is no longer suppressed and final states A* 
include the ground state of the nucleus. The cal- 
culation of the total p spectrum of quarks in qA 
scattering proceeds as follows. The master for- 
mula is 



(PpcPb " (2 



J 



) J At 



(A|Tr{[l - Sa{b)]\A*){A*\[t Siib - c)]}\A) 



(A|Trl|^) 



(25) 



The nuclear S-matrix is a product of the free- 
nucleon S-matrices. The technique of the calcu- 
lation of nuclear matrix elements in the standard 
picture of a nucleus as a dilute uncorrelated gas of 
color-singlet nucleons is found in [Ij, liy and need 
not be repeated here. We simply cite the results. 



{A\TT{Sam\A) 



= S[b,aa{xo)], 



(A\TTt\A) 
{A\TT{Sa{b)Si{b-c)}\A) 
{A\Tt1\A) 



= S[b,aaa{xo,c) 



(26) 



Now we apply H23|l and notice that 

S[b, <Jaa,o{^o)] = S^[b, cr„(xo)], (27) 



which leads to 
(Pp<Pb 



{l-S[b,aa{x,)]}H^^\p) 

+ (l)a-a{b,Xo,p). (28) 



The first term describes the elastic scattering, 
here the delta-function (5(^^(p) is a simple work- 
ing approximation to the nuclear diffraction peak. 
The second term in the p spectrum (I28II de- 
scribes the Born (B) approximation for the inclu- 
sive quasielastic scattering a A —>■ aX, 



dcTQelM 



i{b,XQ,p), 



cPpcPb 
where (j)aa{b,xo,p) is defined in terms of 



S[b,aaaixo,r)] 



C 

S[b, 7^cr(a;o,r)]. 



(29) 



(30) 



Note, that the color representation dependence of 
the dipole cross section (|12|) entails a variety of 
the so-defined collective nuclear glue ^aa{b, xq, k) 
— the nuclear glue is a density matrix in color 
space. The major point is that to the considered 
Born approximation quasielastic scattering is lin- 
ear A;_L-factorizable in terms of the collective glue 

4>aa{b,Xo,p). 

Upon the integration of nuclear Born cross sec- 
tions over all transverse momenta 

a,i = j d?h{l-S[b,aa{x^)]Y, 

(Trn = <JQel^ I d^b{l-S^[b,aa{xo)]}, (31) 
Crtot = CTel + (Jin = 2 d'^b{l - S[b, Cra(xo)]}, 

(32) 

which are the standard Glauber formulas for the 
the projectile with af^ = cra(xo) j4ll |. 



III. MASTER FORMULAS FOR THE 

REAL AND VIRTUAL RADIATIVE 

CORRECTIONS TO INCLUSIVE SINGLE 

LEADING JET PRODUCTION 



A. Real production a 



ag 



The leading order pQCD radiative correction 
to the spectrum of leading particles in aN{A) 
scattering, a — q,g, comes from the radiation of 
gluons, aN{A) -^ agN*{A*). To the lowest or- 
der in pQCD, the generic partonic subprocess is 
agt — > be. From the laboratory frame standpoint, 
it can be viewed as an excitation of the perturba- 
tive \bc) Fock state of the physical projectile \a) 
by one-gluon exchange with the target nucleon. In 
the case of a nuclear target one has to deal with 
multiple gluon exchanges which are enhanced by 
a large thickness of the target nucleus. Here the 
frozen impact parameter approximation holds, if 
the coherence length Ic is larger than the diameter 
of the nucleus 2Ra, 



L 



where 



2Ea 



1 



iQ*)^+Ml xrriN 



> 2Ra , (33) 



^ Zb Zc 



(34) 



is the transverse mass squared of the be state, Pf, ^ 
and Zh.c are the transverse momenta and fractions 




FIG. 4: The color dipole structure (Ihs) of the generic 
3-parton state and (rhs) of the 3-parton state entering 
the color-dipole description of fragmentation a ^ be 
with formation of the be dipole of size p. 



of the the incident parton's momentum carried by 
the outgoing partons (zb + Zc — 1). Notice an 
equivalence between the coherency condition (|33|l 
and the parton fusion condition H19|) . The virtu- 
ality of the incident parton a equals (Q*)^ = k^^, 
where ka is the transverse momentum of the par- 
ton a in the incident proton. 

The target frame rapidity structure of the con- 
sidered a —^ be excitation is shown in fig. 1. The 
beam parton has a rapidity rja >j ija = log 1/xa, 
the final state partons too have rapidities rjb^c > 
r]A- In this paper we focus on the lowest order 
radiative corrections without production of more 
secondary partons in the rapidity span between rja 
(r?,/) and 77c iVg)- 

To the lowest order in the perturbative tran- 
sition a —f be the Fock state expansion for the 
physical state \aphys) reads 



^phys / 



*(z6,p)|6c)o 



(35) 



where \I'(zb,p) is the probability amplitude to 
find the be system with the separation p in the 
two-dimensional impact parameter space, the sub- 
script "0" refers to bare partons. The explicit 
dependence of the wave functions on the virtu- 
ality (Q*)^ and their relation to the parton split- 
ting functions is found in |17J . Here Za is the 
wave function renormalization, the perturbative 
coupling of the a -^ be transition is reabsorbed 
into the lightcone wave function '^{zb,p). The 
normalization (completeness) condition reads 



\^{zb,p)f(fpdzb 



(36) 



If 6q = 6 is the impact parameter of the projec- 
tile a, then 

bb = b + ZcP, bc^b- Zbp, (37) 

see Fig. 0] Below we cite all the spectra in the 
a-target collision frame. We shall speak of the 
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FIG. 5: Typical contribution to the excitation amplitude for aA — > bcX, with multiple color excitations of the 
nucleus. The amplitude receives contributions from processes that involve interactions with the nucleus after 
and before the virtual decay a ^ be which interfere destructively. 



produced parton 6 - or jet originating from this 
parton - as the leading one if it carries a large 
fraction of the beam lightcone momentum, Z{, ~ 1. 
The Feynman variable xp spectra in pN,pA col- 
lisions are obtained from the Zb spectra in qN, qA 
collisions by an obvious A:^-factorization convolu- 
tion with the beam quark densities, see below Eq. 

By the conservation of impact parameters, the 
action of the S'-matrix on \aphys) takes a simple 
form 



^ l^phys) 



J 



Saib)y/Za\a)o 

Sbibb)Sc{bc)^izb,p)\bc)o 
Sa{b)\aphys) 



[Sb{bb)S,{b,) - Sa{b)]^{zb,p)\bc). 

(38) 



In the last line we explicitly decomposed the fi- 
nal state into the elastically scattered \aphys) and 
the excited state \bc)o. The two terms in the lat- 
ter describe a scattering on the target of the be 
system formed way in front of the target and the 
transition a —^ be after the interaction of the state 
|a)o with the target, as illustrated in Fig. \S\ The 
contribution from transitions a —^ be inside the 
target nucleus vanishes in the high-energy limit of 

X <XA- 



The probability amplitude for the two-parton spectrum is given by the Fourier transform 

/ (fbbd^bcexp[-i{p^bb+p^bc)][Sb{bb)Sc{bc)-Sa{b)]-^{zb,p). 

The differential cross section is proportional to the modulus squared of H39() 
da{a* -> 6(pJc(pJ) 



(39) 



dzd^p d^p ' ("2^ j '^^^bd^^c'^^^'bd^^'c exp[iPb(b6 - b[) + ip^{b, - b'J] 

x^bcizb, bb - bc)*L(^6, K - bc){s^tL(^fc' ^- bb, be) + Si'j(6', 6) - S^^Jb, b^ b',) - S^^ib', bb, 6.)}. 

(40) 



where, see Fig.|Slfor an illustration, 

Si^j(6',6) = Sl(6')Sa(b), 
Sgiib\bb,b,) = Si{b')Sb{bb)Sc{b.), 



Smb,blb'J - S^(6;,)St(6^)S„(6), 



bca 



Sfl^{blb'„bb,bc) = Sl{bi)Sl{b'JS,{b,)Sb{bb) . 

(41) 



Here the proper averaging over polarizations 
and color states of the initial parton and the sum- 



11 






— c 

H 
C 

b 






K-< 



bccb aa 



q(3)_ 
^acb 



<: 



q(_3) 
^abc 



FIG. 6: The S-matrix structure of the two-body density matrix for excitation a — > be. 



mation over polarizations and color states of final 
state partons are understood. Next we approxi- 
mate the heavy nucleus by a dilute gas of colorless 
nucleons and follow Glauber's technique of sum- 
ming over the final states of the target nucleus ap- 
plying the closure relation |4l|. Then the nuclear 
matrix element of products of single-parton S- 
matrices would boil down to the Glauber-Gribov 
formula ^■. (A|S("^|A) = S[b;E(")], where S(") 
is the relevant n-parton cross section operator, for 
the details see Refs. [13, El Hi. All the multi- 
parton states - aa, abc, bca, bccb - are in an over- 
all color-singlet state and the excitation cross sec- 
tions are infrared-finite despite the nonvanishing 
net color charge of the projectile parton a p!7ll37l|. 
Integration over the transverse momentum p^ 
of the jet c gives 

K = b'^, b-b' = zt{p-p'), 
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p p', b'-b, = Zbp'. (42) 



Then the unitarity relation 

SUbc)Sc{bc) = 1 
leads to a fundamental simplification: 



(43) 



^htrh(^b^^'c,^b,bc) 



- Sl{b',)Sl{b,)S,{b,)St{b,) 

= Slib',)St{h) 

= si^Ublb,), (44) 



i.e., the effect of interactions of the spectator par- 
ton c vanishes upon the summation over all its 
color states and integration over all its transverse 
momenta (|33|, see also a discussion in [131) the 
only trace of the observed parton b having been 
produced in the fragmentation a — > 5c is in the 
density matrix '^(zfc, bf, — 6c)^*(z6, b^ — b'^) which 
defines the transverse momentum distribution of 



the parton b in the beam parton a and the parti- 
tion Zf, , Zc = (1 — Zh), of the longitudinal momen- 
tum between the final state partons. The applica- 
tions of the master formula (|4()|l to the calculation 
of the real emission spectra for the free-nucleon 
and nuclear targets are found in |l7| . 



B. Derivation of the virtual radiative 
correction from the unitarity relation 

The relevance of the unitarity has already been 
explained in the Introduction. It is especially ob- 
vious in qN{A) scattering. The opening of the 
new channel q —^ qg renormalizes, via unitarity, 
the Born cross section of the lower order quasielas- 
tic scattering qN{A) -^ q'N*{A*). The VRC 
is precisely that unitarity driven renormalization. 
The real emission and the virtual correction must 
so conspire, that the total multiplicity of leading 
quarks in the final state is exactly unity. We show 
here two derivations which have already been used 
in the related works [33, [sl, [S^l ■ 

The first approach is a direct perturbative solu- 
tion of the unitarity relation (J37l|. for an extensive 
use in the theory of the LPM effect see [if). We 
define the radiatively corrected S-matrix of elastic 
scattering, 

Sei{b,x) ^Sa{b,xo) + SSa{b,x), (45) 

and the excitation (real emission) operator 

Sex{b,x) ^Sbibb,xo)Sc{bc,XQ) ~Sa{b,xo). (46) 

Then the unitarity relation would read 

S{b,x)S^{b,x)=Sei{b,x)Sl{b,x) 

+ fvbcn,^bcSeAb,x)SUb,x) = 1, (47) 
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where we explicitly separated the contributions the interval [a;/a;o: !]■ We recall that ^J^^bc ex a^ 

from the elastic and real emission channels, Pbc is and the real emission contribution to H47|l must 

the short hand notation for the integral over the be regarded as the pQCD perturbation. Making 

be dipole parameters, Dbc — (Ppdzc- If the bound- use of the explicit form of Sei(b, x), to the linear 

ary condition for the small- a; evolution is defined order in dSa{b,x) we find 
at xq, then the integration over dzc extends over 

I 



SSa{b, x)Si{b, xo) + Saib, xo)SSi{b, x) = 2?bc*6c*bc[Sfc(6b, a;o)Sc(6c, a;o)Sl(6, xq) ~ 1] 

+ / Vbc%,-^bc[Sa{ba, xo)Sl{b,, xo)Sl{bb, Xo) - 1]. (48) 

Evidently, this equation can be split into 

6Saib, x)Siib, xo) = f Vi,cn,-^bc[Sbibb, xo)Sc(bc, xo)Si{b, xo) - 1] (49) 

and its complex conjugate. Multiplication of H49|) by Sq(6, xo) from the right and the application of the 
unitarity relation (|10|l gives the solution 

6Sa{b, x)= Pfcc*hc*bc[Sb(bh, xo)Sc{bc, xq) - Sa{b, Xo)] (50) 

The second approach is a direct calculation 

Sei{b,x) = {aphys\S\aphys) ^ ZaSaib,xo) + Vbc'^l^^bcSb{bb,xo)Scibc,xo) 

= Sa{b,xo)+ Vbc-^*b^^bc[Sbibb,xo)Sc{bc,xo)-Sa{b,xo)], (51) 

where in the last line we applied the completeness condition H^fifl . It gives precisely the same result 
for SSa{b, x) as the unitarity condition. The technique (|51|l was applied earlier to the derivation of the 
color dipole BFKL equation ^, ^J . 

IV. THE CALCULATION OF THE VIRTUAL RADIATIVE CORRECTIONS TO THE 

LEADING PARTON SPECTRUM 

A. Leading partons off the free-nucleon target 

With the allowance for the radiative correction, the quasielastic scattering matrix 1 — Sa(b, xo) in the 
master formula l|16(l must be replaced by 1 — Sg; (b, x) and 

[1-Sel{b,x)][l-Sl{b',x)] = [l-Sa{b,Xo)][l-Si{b',Xo)] 

- SSaib, x)[l - Sl(b', xo)] - [1 - S,(6, xo)]6SUb\ x) (52) 

Here first term leads to the Born spectrum. Upon using of the result H5()(l . the virtual correction, 
j^^VRC ^ -j-p |.]-^g Born spectrum takes the form (hereafter we suppress the calculation of the free-nucleon 
and nuclear matrix elements) 

da^^'^(x) f d'^bd^c f 

rf2 " "^ '/ (27r)2 °^P("*P • ^) / 2?fcc*fcc*b4S6(bfc, xo)Sc{bc, xo) + 

Sq(6, xo)S'l{b - c, Xo) - Sa{b, xo) ~ Sbih, xo)Sc{bc, xo)Sl{b, xo)] 
d^c 



(27r) 



exp(-ip • C) / Vbc%c^bc[crabcixo) + (7a{xo) ~ (Tbcixo) - CTaaixo)] (53) 
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The general technique of the derivation of the multiparton cross sections is found in (Tfl, lO . ll4L l3a | , 

here we only cite the results: 



(Tbc{x) 



-cr(x,p) + 7-^^0-0(2;), 



2Cf 



2Cf 



cabcla;) = ^r;^ cr(x, bf, - be) H ^r^^ (^[x,b - be) ^ — ; a{x,b - bt) 



2Cf 
2CV 



2Cf 



2Cf 



I \ , ^a + L/C ^ L/b / , N , ^a + L/ft — Gc , , 

cr(x, p) H — ; cr^a:, c + Z{,pj H -— a[x, c - ZcP). 



2Cf 



2Cf 



(54) 



The integrand of (|53|l takes the form 

(yabc{x) + (Ta{x) - abc{x) - (7aa{x) = 

-[{Ca + Cc- Cb)cr(x, c + Zbp) + {Ca + Cfc - Cc)CT(a;, c - ZcP) - 2Ca(T(a;, c)] 



2Cf 
1 

2C^ 



— / (fKf{x, k) exp(JK • c)[2Cq - {Ca + Cc- Cb) CXp{iZbK ■ p) - [Ca + Cb- Cc) cxp(-iZcK • p)]. 

(55) 
Now we specify the phase space limits. Let x — p'^ /W^jy and the boundary condition be specified at 



xq. Then the Zc integration goes over [zmtn, z^nax], where z^nax = 1 - Zmin, and ZmaxZTmn 

The resulting VRC to quasielastic scattering can be cast in two forms. In the color dipole representation 



.xoW^^ 



p\ 



da^^'-{x,p) 
<Pp 



1 



2Cf (27r)2 



fixo,p) 



dZc / d^p^le{Za,p)^bciZc,p) 



2Ca - {Ca + Cc- Cb) exp{izbP ■ p) - {Ca + Cb~ Cc) exp(-iZcP • p) 

dcTQclM Ca 



^ / dZc / dV*6c(^c,p)*6c(^c,p) 



d^p Cf (27 

|[1 - exp(izf,p • p)] + [1 - exp(-izcp • p)] + (— ^ - l)[exp(-iZcP • p) - exp{izbp ■ p)]j 

(56) 



has the expected form of the Born cross section times the combination of form factors of the be Fock 
state. This combination of form factors vanishes at p = 0. In the final result we made an explicit use 
of Cb = Ca and Cc — Ca which is the case in reactions of the practical interest: a -^ ag. 

In the subsequent analysis of the small- cc evolution properties of the leading particle spectrum it would 
be more convenient to work in the /c^-factorization representation 



da^«^(x,p) 1 

= f{xo,P) 



1 



d^p 



dzc d K 



2Cf 2(27r)2 

+ {2Ca~CA)\'^{Zc,K)--i'{Zc,K + ZcP)\'' 



Ca\'^{zc,k) -^{Zc,K- Zbp)\ 



1 1 /"^ 

/(a^cP):^;:^ • :77T7^ / dzc2asPca{zc) 
2Cf 2{2tt)^ 7^/^^ 

/ d^K[CAK{K, K - Zbp) + {2Ca - Ca)K{k, K + Zcp)] 



(57) 



I 

The averaging over the initial, and summing is encoded in [^^(zc, k) — 'I'(2:c, /« —p)P, for which 
over the final, colors and polarizations of partons 
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we use the explicit form [T3 

|*(zc,k)-*(zc,k- A)p = 

= 2asPca{Zc) 



K 



K^ 



£2 (/«-A)2+£2 



2asPca{Zc)K{K,K-A) 



where 



ZbZcQl + Zbml + Zc 



% 1 



(58) 



(59) 



(32 is the virtuaUty of the projectile parton de- 
fined by its transverse momentum in the incident 
proton, (52 = /g^^ Q^Yifi Pcaizc) is the real-emission 
part of the familiar splitting function. For inci- 
dent quarks, a — q, 



^ga\Zg) ^a 



1+(1-^,)^ 



(60) 



for incident gluons one would use the splitting 
function Pgg{zg), which also will be relevant to 
the color octet incident quarks. We emphasize in 
passing that in transitions a — > 6c the helicity of 
partons mixes with the orbital angular momen- 
tum, but that has no bearing on the description 
of the polarization summed final states because 
the s-channel helicity of partons is conserved in 
the scattering process. 

In the calculation of the total radiative correc- 
tion, the above VRC must be combined with the 
(nonlinear) k factorization results for the real 
emission [l3|. Making use of the representation 
lf55|) . the spectrum of color-triplet leading quarks, 
b — a = q, from the real emission a ^> bg can be 
cast in the form (p = Pf,) 



da^'^(^{x,p) 



dzhCpp 



a^bg 



2asPga{Zg) 



2Cf 2(27r)2 

/ (fKf{xa,K.)\CA[K{p,p-K) 

K{p,p- Zbn)] 

{2CF-CA)K{p~K,p-ZbK)y (61) 



The corresponding result for the color-octet inci- 
dent partons is 



da^^^{x,p) 



dzbd^p 
1 1 



a^bg 



■2asPga{Zg 



2Cf 2(27r)2 

X Ca d'^K,f{xQ,K)lK{p,p-ZgK) 

+ K{p-K,p-ZgK) + K{p,p-K)y {62) 

Once the Rnite-Zg VRC and REC for quarks 
are known, we can proceed to inclusive jets in 
collisions of protons. Here a crucial point is a 
cancellation of effects of interactions of spectator 
partons of the beam proton in the inclusive single- 
particle cross section summed over all excitations 
of the beam and target [13, 123 ■ Let qa{xN,ka) 
be an unintegrated density of quarks of fiavor a, 
carrying a fraction Xa of the proton's lightcone 
momentum and having the transverse momentum 
ka- Then the inclusive spectrum of leading jets 
with the Feynman variable xp and the transverse 
momentum p, produced by the interacting quark 
a will be given by the standard convolution, cf. 
Eq. ®, 



daa{xF,p) 
dxpd'^p 



/ dXa d'^ka dZbd^p'5{xF - ZbXa)5{p-p' ~ Zbka) 

JxF J Jo 



da^Jx,p') da^"^{x,p') 



d^p' 



d^p' 



<5(1 - zb) + 



da^^^ix,p') 
dzbd^p' 



(63) 



where XaxWppf = p^ / ZbZc and WpN is the proton-nuclcon cms energy. 



B. Color-triplet leading quarks off the nuclear target 

In the case of heavy nuclei one considers cross sections per unit area in the impact parameter plane. 
Starting with Eq. H53|) . we calculate first the nuclear matrix element of the multiparton S-matrices. 
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The resulting master formula for the VRC takes the form 



(PpcPb 



-S[b;aabc{xo)]} (64) 



For color-triplet incident quarks Ca — Cb ^ Cp- Making use of (|54|l . and treating the nucleus as the 
dilute uncorrelated gas of colorless nucleons, we readily find the product representation for multiparton 
S-matrices: 

S[b; (Tbc{xo)] + S[6; aaa{xo)\ - S[b; o-a(a;o)] - S[b; (Jabcixo)] 

= {S[b;^^fT(a;o,p)]S[5;-cro(a;o)] - S[b; -ctoK)]} 



Ca 



2Cf ~ Ca 



|S[b; a{xo,c)] - S[b; ^^cr(a:;o, p)]S[6; -^(j{xo,c + Zbp)]S[b; — ^ — -(j{xq,c - Zgp)]^ 



2Cf 



2Cf 



2Cf 



(65) 



The Fourier transform of the first group is of terms, which is c-independent, will be oc 5^ \p). We 
evaluate the second group to the leading order of the large- A'c perturbation theory, when Ca = 2Cf 
and S[b; jrT '^(^Oi c — Zgp)] = 1. In this approximation it admits a simple Fourier representation 



S[b; (t{xo, c)] - S[6; -— -cr(xo, p)]S[b; -— -cr(a;o, c + Zbp)] 



2Cf 



2Cf 



/ d^K(f>{b, xo, /«) exp(iK ■ c) x <1 — S[6; cr(p)] exp{izbK ■ p)f. 



(66) 



The VRC splits naturally into two terms 



d^pd^b 



= 2(5(2) (p) / dz,d^pnc{Zc,p)-^bc{Zc,p) 



X {1 - S[6; cr(a;o, p)]} • S[b; -^0(2^0)] • {1 - S[b; -cro(a;o)]} 



2(/)(b,a;o,p) 



dzcd^ p%^{zc, p)^bc{zc, p){l - S[b; a{xo,p)] exp(izfcp ■ p)}. 



(67) 



Evidently, the term oc (J'^-* (p) in the r.h.s. of Now, the VRC to the profile function of pure elas- 



(jB7|) is the VRC to the elastic scattering. Indeed, 
the VRC SSei{b) to the elastic S-matrix equals 

SSelib) = {A\SSa{b)\A) 

= dZcd^p^l^{Zc,p)^bc{Zc,p) 

x{S[b;abc]-S[b;<Ja]}. (68) 

With the large- iVc cross section (|54ll . this gives 

5Sei{b) = -S[b; -(To{xn)] 



dzcd^p%^{z^, p)^bc{zc, p) 
X {l-S[b;a(xo,p)]}. (69) 



tic scattering can be evaluated from the expansion 



{l-S[b,-a^{x^)]-5S,i{b)Y 

= {l-S[b,iao(xo)]r 
2,5Sez(b){l-S[6,i(7o(xo)]}. (70) 



In conjunction with (|69|l one obtains precisely the 
integrand of the term oc 5^'^\p) in the r.h.s. of 
(|67|l which accomplishes the proof. 

Now we isolate the nonlinear k±^ factorization 
for the virtual correction to quasielastic scattering 



of color-triplet quarks: 

d?p(Pb 
- 2(f>{b,xo,p) / dzc(f p^l^izc, p)-^bc{zc, p) 

X {l-S[b;a{xo,p)]exp{izbp- p)} 

= -0(b, xcp) / dZc2asPgq{Zc) / d^Kl 

X lS[b\(To{Xf))]K{K.i,Ki- Zbp) 

d'^K2(t>{b, Xo,K2)K{ki,Ki + K2 



zbP)Y 



(71) 

Here we made use of the Fourier representation 
(|23|) for S[6; cr(p)] and the sum rule H24|) . In view 
of Eq. (|29|l . and in close similarity to the free- 
nucleon case, the VRC has the expected form of 
the Born cross section times a p dependent form 
factor. The difference from the free-nucleon case 
is that this from factor does not vanish at p = 0. 
The above VRC must be combined with REC 
derived in Ref. [I3| to the same leading order of 
the large- -/Vj, perturbation theory {a — b ~ q): 



{2^fda^^^(x,p) 



dzbd'^pd'^b 

/ d^Ki(p{b, Xo, Ki 



2asPgq{zg) 



q^gg 



S[b; aoixo)][K{p,p- k) + K{p,p - zi,k)] 
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C. Color-octet leading quarks and leading 
gluons off nuclei 

The case of color-octet incident partons gives 
still more insight into the leading parton produc- 
tion. To the Born approximation, the spectrum 
of color-octet leading partons (|29|l is described 
by the glue 4'gg{b,XQ,n) defined through the S- 
matrix for the octet-octet color dipoles \u\ : 



$ 



gg{b,XO,p) = 



(Pr ^r, Ca , ., 



{2nf ' ' Cf 
xexp{—ip-r). (73) 



The issue is how this result changes if one in- 
cludes the VRC and REC for gluons radiated in 
the regime of nuclear coherency in the rapidity 
span between the jet and the target nucleus. 

A tricky point is that in our derivation we shall 
encounter not the collective glue (ff!^ but still an- 
other component of the color-space density matrix 
for their collective glue. Specifically, in the case of 
incident color-octet partons Ca — Ct ^ Cc = Ca 
and all the multiparton dipole cross sections (|54|l 
are superpositions of 



Ca 



1 



'^g{x,r) = ^^cr{x,r):^ -agg{x,r), 

f^gflix) = ^^^o(a;) = -j^(^ggfl{x). (74) 

Correspondingly, the multiparton nuclear S- 
matrices will be products of S[b;Ug(x,c)\, and 
there emerges a collective nuclear glue $g(b, x, k) 
such that 



^ j d'K2mx,.i^2)K{p^K,.^p~Zbn2))^. S[b;ag{x,c)]^ J d'K<i>,ib,x,K)eM^>^-c) 



(75) 



(72) 

For a practical calculation of the LPM effect and 
nuclear quenching of leading quark jets off nuclei 
one would use an obvious generalization of the 
convolution l(^ to nuclear targets. 



J 



In contrast to ^{b,xo,p) and ^gg{b,xo,p), the 
so-defined auxiliary <I>g(6, xq, k) is a component 
of the color density matrix for the collective nu- 
clear glue which is not directly measurable in the 
quasielastic scattering of partons. 



What would enter the master formula for color-octet leading partons is 



S[b; abcixo)] + S[b; aaa{xo)] - S[b;aaixQ)] - S[b; aabc{xo)] 
= S[b; (Jg{xQ, p)]S[b; agfiixo)] - S[b; ag^o{xo)] + S[b; 2ag{xQ,c)] 
-S[b; <7g{xo, p)]S[b; <7g{xo,c + Zbp)]5[b; ag{xo, c - Zgp)] 
= -S[b; (Tg,o(a;o)]{l - S[b; crg,o]}-{l - S[6; crg(xo, p)]} 



-S[b; 0-3,0(2^0)] / d^Kexp{iK ■ c)(/)g(b, a;o/«)|2 - S[b; crg(xo, p)] • [exp(izf,K • p) + exp(-iZgK • p)]\ 
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+ / (fKi(fK2<l)g{b,Xa,Ki)(l)g{b,Xo,K2)ex.p{i{Ki+K2)-c) 

xjl - S[b;(Tg(a;o,p)] exp[ip(zf,Ki -iZgK2)]\. 



(76) 



Again we made an extensive use of the dilute gas treatment of the nucleus. Precisely as it was the case 
with the color-triplet quarks, the c independent term in ifT^ gives rise to the VRC to the pure elastic 
scattering: 



(Pp(Pb 



= (5^2) (p) . 2S[6; cTg^oixo)] ■ {1 - S[b; crgfi{x^)]} / dz^d^ p%^{z,,, p)^bc{zc, p){l - S[6; (jg{xo,p)]}. 

(77) 



The nonlinear k± factorization for the quasielastic scattering of the color-octet quarks is much more 
subtle: 



(27r) 



da^'^^jx^p) 



dzq2asPc 



,a{Zg) / rf^Ki 



X < (t>gib, XQ,p)S'^[b; ag^Q{xo)][K{Ki, Ki + z^p) +K{k.i,Ki + Zgp)] 

- S [6; CTg, o(a;o)] d'^K2l(j>g{b,Xo,p-H2)(f>g{b,Xo,K2)K{Ki,Ki-K2 + Zgp) 
+ (l)g{b,Xo,p)(j)g{b,Xo,H2)[K{Ki,Ki + H2 + Zbp) + K{ki,Ki +K2 + Zgp)]\ 

- / d^K,2d^K3(l)g{b,Xo,p- K2)4'g{b,Xo,K2)(l)gib,Xo,K3)K{Ki,Ki - K2 + K3 + Zgp) >. (78) 

In striking contrast to the case of color-triplet quarks, the VRC to quasielastic scattering of color-octet 
quarks is no longer proportional to the Born approximation. In the calculation of the total spectrum of 
color-octet leading quarks, the VRC H78|) must be combined with the REC derived in Ref . [13 : 



{2Try da'^^^ixp) 



dzbd^pd^b 



a^bg 



= 2asPga{zg) / d^Ki<j)g{b,xn,K.i) 



xl S'^[b;(Tg^o{x())][K{p,p + ZbKi) + K{p + K,p + ZbK) + K{p,p + K)] 
-|-S[b;crg,o(2;o)] / d^K2(j)g{b,XQ,K2) 

x[K(p + Ki,p+ ZbK2) + K{p+ Ki + K2,P+ ZbK2) + K {p + Ki,p + Zb{Ki +K2))] 
+ / d'^K2d'^K3(j)g{b,Xo,K2)(j)g{b,XQ,K,3)K(p + Ki +K3,p + Zb{K2 + /^s)) }■ 



(79) 



r 



D. The s-channel unitarity and the quark 
number sum rule 



Although our derivation of the virtual radiative 
corrections makes manifest use of the s-channel 
unitarity, it is useful to revisit this issue from the 
viewpoint of the leading quark-number sum rule. 
The latter dictates that the leading quark inclu- 



sive cross section must integrate to exactly the to- 
tal quark-target interaction cross section, i.e., the 
sum of the integrated virtual and real radiative 
corrections to the single-quark inclusive cross sec- 
tion must equal the increment of the total quark- 
target cross section caused by the gluon produc- 
tion in the s-channel. 

We start with the integrated real emission cross 
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section. The integration over the transverse mo- 
menta p^ and Pfj in the master formula H40() gives 
be = b'^ and bi, = b[, which entails b = b', 
p — p' . Then, the unitarity condition H43|) 

leads to S^'tl^{b[,b'e,bb,bc) = 1 and S^^j(6',6) = 
1. The three-parton cross section which enters 
S\{b')Sbibb)Scibc) would equal 

O-abc = TZ7^ cr(p) H (j{zbp) 



2Cf 

Ca + Cb — Cc 

2CZ 



2Cf 



a(zcp). 



(80) 



The result for the integrated real emission cross 
section is (we suppress the matrix elements over 
the target nucleon and/or nucleus) 

5a'''''' (a* ^ h{pb)c{pS) = 2 j Vbcnc^bcd^b, 

X[l-Sahc], (81) 

cf. the derivation of open charm production in 
Ref. [13. 

The integration over the transverse momentum 
p in the master formula (|53|l amounts to c = 
and 



5a^^^ = -2 / mL^brVbrd^b. 



bc^bcl^bcl^ 

^[l + Sbc-Sa-S-abc]. (82) 

Then, the total radiative correction to the in- 
clusive single-quark cross section would equal 



5a 



REC+VRC 



2 I %^^bcVbcd^b,[Sa-Sbc] 



= -2 / n^-^bcVbcd''b,5Sa 

= 5atot- (83) 

In the last step we used our result (|66|l . We em- 
phasize that Eq. l|SS|) gives the exact result for 
the increment of the total quark-gluon cross cross 
section for radiation of one gluon, without resort- 
ing to the soft gluon, i.e., LLi approximation. 
The proof of the quark-number sum rule is ac- 
complished. 



E. A mini-summary on the LPM effect for 
leading jets and particle 

The nuclear modification of the radiation 
quenching (stopping) of leading jets produced at 
large Feynman variable and fixed transverse mo- 
mentum - the p-dependent LPM effect for lead- 
ing jets - is our principal new result. The non- 
linear k\ factorization of Sees. IV. B and IV. C 



gives the leading particle spectra in the form of 
explicit nonlinear quadratures in terms of the col- 
lective nuclear glue. Such quadratures for both 
real and virtual pQCD radiative corrections are 
not found in the previous works by other groups 
on forward jet production [l^ [l^ [23] ■ Still an- 
other novelty of our approach is a derivation of 
the virtual pQCD radiative corrections from the s- 
channel unitarity. Our approach to the solution of 
the s-channel unitarity holds for both free nucleon 
and nuclear targets. The emerging interplay of the 
real and virtual radiative corrections to the lead- 
ing jet spectra is precisely the same as in the evo- 
lution equation for the total cross section which 
property guarantees a manifest fulfillment of the 
leading parton number sum rule. For color octet 
(adjoint representation) leading quarks and gluons 
our nonlinear k±^ factorization holds for arbitrary 
A^cj only for color-triplet (fundamental represen- 
tation) leading quarks one would invoke large- iVc 
perturbation theory for the formulation of closed- 
form quadratures. 



V. LINEAR fcx FACTORIZATION FOR 

THE LEADING PARTON SPECTRA TO 

LLi? 



A. The master formula for LL^ evolution of 

the color-dipole S-matrix and the collective 

nuclear glue beyond the Born approximation 

For further insight into the properties of lead- 
ing parton spectra, in this section we consider the 
LPM effect-deconvoluted cross sections, i.e., the 
spectra evaluated in the LL-, i.e., soft gluon, ap- 
proximation neglecting the radiative energy loss. 
The relatively simple case of the free-nucleon tar- 
get will be considered in the next subsection, only 
with the nuclear targets one encounters nontriv- 
ial conceptual issues. Here, motivated by the lin- 
ear fc_L-factorization H29(l observed to the Born ap- 
proximation, we would like to compare the LL^ 
evolution properties of the leading particle spec- 
tra and of the collective nuclear glue. We recall 
that in Ref. [H HI HI we defined the collec- 
tive nuclear glue in terms of the amplitude for 
coherent diffractive dijets. At arbitrarily small 
X this observable is calculable through the color- 
dipole S-matrix Saa{x,ba.,ba)- At a boundary 
xq = xa, i.e., at the onset of coherent nuclear 
effects, Saa{xo,ba,ba,) is given by Eq. if^ which 
defines the Born approximation for the collective 
glue (|23|l . We would like to maintain the def- 
inition l|23|) for arbitrary values of x and need 
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the LLi evolution properties of the color dipole 
S-matrix. For soft gluons, Zg -^ 1, there are 
several important simplifications. First, the im- 
pact parameters of partons which radiate soft glu- 
ons, do not change, see Eq. l|H7|) . The resulting 



master formula for the small- a; renormalization of 
the color-dipole S-matrix for the perturbative aag 
Fock state in the a a-dipole is an obvious general- 
ization of Eq. l|CT|l 



g-gi poie 

Mm- 



SSaaix,ba,ba) ^ / dZg (f p'i'l^g'^ aag[Saixo, ba)Sl{xo, ba)Sg{xo, bg) - SaaixQ, ba, ba)]. (84) 

Jx/xq J 



Second, in view of Eq. 159|l . the parameter e^ in agation radius of perturbative gluons 39,^44]. The 
the ag wave function would not depend on the corresponding small- Zg splitting function equals 
virtuality of partons a, a in the aa color dipole. 
Consequently, the distribution of gluons in the aa 
dipole will be given by 38, 39] 



or' 

PgaiZg)^^^. (87) 



Kag'^a-ag = |*ag(^<„ p) ^ ^agiZg, P + r)f, (85) 

where r = ba — ba, p = bg — ba- In the momentum 
space 



"^agiZg^k) CX -^- -, 



(86) 



To the LLj we have [zmin,Zmax] = [a;/a;o,l] and 
the Zg integration would give the familiar LLi 
factor 



dzg XQ 

- log—. 



where we indicated the (optional) infrared regular- 
ization parameter ijg which models the finite prop- 



x/xo ^9 



J 



B. Free-nucleon target: BFKL evolution for the leading parton spectrum 

We start with the free-nucleon target. We first recall the connection between the LL^ evolution 
equation for color dipole cross section [33, |33 and the BFKL equation 0, Q for the unintegrated glue 
of the nucleon. 

In terms of the dipole cross sections, Eq. l|5H) amounts to the evolution correction |33, 13^ 



6aaa{x,r) 



/ dZg j d^p\^ag{Zg,p) -■^ag{Zg,p + r)\'^{cFaag{xQ) - CFaaixo)) 

^ JdZgJ d^pl'^agiZg, p) - ^ag{Zg. P + r)\'' 



X [o-(a;o,p) +CT(xo,p + r) -o-(a;o,r)], 
and to the evolution of the unintegrated glue 6f{x,p) 
d^r 



Ca 



.„ .^Saaa{x,r) exp{-ip ■ r) = Saaa.o{x)S^'^Hp) - -^Sf{x,p). 

(ZTTJ Up 



(89) 



(90) 



In the transition from the color-dipole to the momentum representation we follow the familiar route. 
First, we make use of 



aix,p)+aix,p + r)-.ix,r)^ld^nfix,.)il^oM^^-pm + cM^'^-r)]. 



(91) 
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Second, the symmetry with respect to p <-^ p + r allows a simplifying substitution of \^ag{zg, p) — 
^ag{zg,p + r)\'^ by 2**g(zg,p)[*ag(zg,p) - -iiagizg^p + r)] which has a Fourier representation 

"^IgiZg, p)[^ag{Zg, p) - "ilagiZg, p + r)] = / (f Ki(f K,2^* {Zg , ^2)^ {Zg , Ki) exp[i{K,i - K2) ■ p] 

x[l-exp(iKi -r)]. (92) 

Then the Fourier transform (|90|l gives 



^a c J., ■, Ca 1 
—Sf{x,p) = 



/ dZg / (pKCpKitpK2f{xo,K)'^*{Zg,K2)'i'{Zg,Ki) 



Cf (27r)2 

X [6(^\K^^n2)-5^^\K + Ki-K2)] 

X [J(2) (ki -p)+ 5(2) (k, + Ki - p) - ,5(2) (k - p)] 



Ca 1 



dzg2asPga(zg 



Cf (2^)2 
/.=»{2A-(p.p-.)/,«,.,-/(..,p,A-K»-p,} 



(93) 



Upon the use of the soft-gluon approximation (p?7jl , the differential form of Eq. IplH) ! boils down to 
precisely the BFKL equation for the unintegrated gluon density of the target nucleon [^ Q : 



df{x,p ) 
Slogi 



^ = ICo f d^K[2K{p,p- K)f{x,K) ~ f{x,p)K{K,K^p)] ^ (^ICbfkl® f){x,p) (94) 



r 



where 



/Co 



2^2' 



(95) 



The manifest dependence on the Casimir Ca of 
projectile partons in the Ihs of Eq. (|93|l dipole 
cancels against C'a in the a —> ag splitting func- 
tion. Similarly, we obtain 



dSaaa,Q{x) 
aiogi 



Ca_ 

Cf 

Ca_ 

Cf 



/Co / d^qd^KK{q,q + K)f{x,K) 

(96) 



a^Jfi-,-) 



91og- 



Note the consistency between equations (|96ll and 

Now we proceed to the LL^ approximation for 
quasielastic scattering. In the REC the energy 
loss by the quasielastically scattered parton is ne- 
glected, i.e., one must put Zb = 1 where legitimate 
and integrate over Zg with the soft-gluon approxi- 
mation 1)87(1 . The contribution to the integrand of 
l|CT)) from the term oc {2Cf — Ca)K{p — k.,p—zi,k) 
vanishes. Similarly, the contribution to the VRC 
from the term ex {2Ca — CA)K{K, k+ZcP) vanishes 
at Zc = Zg ^ 1. Combining together the virtual 
and real corrections and the Born cross section, 



we find the z-integrated inclusive p spectrum of 
leading partons 



daQei{x,p) _ Ca 



^lfixo,p) + JColog^Jd^^ 



d^p 2Cf 

X [2K{p,p - K)f{xo,K) - f{xo,p)K{K, K - p)] 



Ca 

2Cf 

Ca 

2Cf 
dcrQei{xo,p) 
d^p 



{ fixo,p) +log — (^ICBFKL®fj{^0,P) 

fix,p) 



log — (iCbfkl (8) -1^) ixo,p)- 

(97) 



We conclude, that to the LL^ approximation, 
the inclusive spectrum of leading partons maps 
the BFKL-evolved unintegrated glue in the tar- 
get nucleon, as it was announced in the Introduc- 
tion, Eq. ^. This spectrum by itself satisfies the 
linear BFKL evolution! It is proportional to the 
Casimir of the projectile parton Ca, i.e., satisfies 
the expected Regge factorization. Of course, these 
Regge factorization properties are precisely what 



one would have expected from the dominance of 
multiregge production processes to LLj lalT'l- In 
our approach, the VRC component of the BFKL 
kernel - the origin of the reggeization of gluons 
- derives from the manifest imposition of the s- 
channel unitarity; such a derivation is of obvious 
relevance to the case of nuclear targets with severe 
unitarity constraints. 



C. The first LL^ iteration for nuclear glue: 

triplet-antitriplet qq dipoles and color-triplet 

leading quarks 

The principal reasoning behind our definition 
|2J, I23 of the collective nuclear glue - the linear 
fc_L-factorization for the amplitude of excitation of 
coherent diffractive dijets - remains valid beyond 



J 
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the nuclear Born approximation at a; « xa, see 
below Sec. V.E. Consequently, we stick to (P5|l 
as the all- a; definition of the x-dependent collec- 
tive nuclear unintegrated glue $(b, x, p). We need 
to specify the LL^ evolution properties of this S- 
matrix. As we saw already at the Born approx- 
imation, the collective nuclear glue is a density 
matrix in the color space. We need to investigate 
how the LLi evolution depends on the color mul- 
tiplet the partons of the color-dipole belong to. 

We start with the first iteration of the LLi 
evolution for the triplet-antitriplet color dipole S- 
matrix to the leading order of the large- A^c per- 
turbation theory in which Ca = 2Cf. The ef- 
fect of the perturbative qqg Fock state in the qq 
color dipole is given by a nuclear generalization of 
Eq. IHl: 



6Sqq{b;x,r) = dzg (fp\'^agizg,p)-'i'ag{zg,p + r)\'^lS[b;aqqg{xo)]-S[b;ao{xQ)]\ 

= dZg (fip\^ag{Zg,p) ~'^ag{Zg,p + r)\^ 

X |S[b;cr(a;o,p)]S[6;cr(a;o,p + r)] -S[6;cr(xo,r)]|, (98) 



where we made use of the large- iVc approximation 
for (jqqg of Eq. H54|) (for the early suggestions to 
consider (|98|) as the closed-form equation for the 
color dipole S-matrix see Refs. 29, 30], the recent 
dispute and further references are found in Refs. 
|23,|43)- The first LL- iteration for the collective 
nuclear glue S4>(b,x,p) is defined by 



S^{b,x,p) 



—-—6Sqq{b, X, r) exp{-ip ■ r) 
SSqq,o{b,x)S^^\p) + S^{b,x,p) 



(99) 

Following the technique exposed in Sec. IV, 
one readily finds the nonlinear fc_L-factorization 
quadrature (the preliminary results were reported 
at conferences |31|) 



d<^[b,x,p)^ / dzg2asPgq(Zg) / — — - 

Jx/xo J \^'^ 



2 



X <^ ^{b,Xo,K,i)^{b,XQ,K2)K{p+Ki,p + K2) 
- I (fKi(fK2^{b,Xo,p)^{b,Xo,fi2) 



xK{ki,Ki +K2) 



(100) 



Upon the application of the expansion (|23() . this 
equation splits into 



SSqq^oib,x) = -S[6;CTo(a;o)]^olog 



Xo 



X d'^K.(fqK{q,q + K)<j){b,xo,n), (101) 

where we notice a nice correspondence to the re- 
sult (PH|l for the free-nucleon target, and 



5(j}{b,x,p) = /Colog — < S[b;cro(a;o)] j d^K 

X 2K{p,p — K)(f){b,XQ,K) — 
(t){b,xo,p)K{K,K -p) 

+ d'^Ki(fK2(f>{b,XQ,K2) 

X K{p+ Ki,p + K2)(f>{b,Xo,Ki) 
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-K{Ki,p+Ki+K2)(l3{b,Xo,p) y 

(102) 
The evolved collective nuclear glue equals 

(j){b,x,p) = (j){b,xo,p) + 5(j){b,x,p). (103) 

Notice, that the result H102|) contains the linear 
component which evolves with the conventional 
BFKL kernel, but this component is suppressed 
by the nuclear attenuation factor S[6; cro(2;o)] 
(see, however, below Sec. V.G). If we identify 

I 



the collective glue (j){b,xo,p) with the exchange 
by a nuclear pomeron, then the absorption-non- 
suppressed component of l|102(l can be viewed as 
a fusion of two nuclear pomcrons. 

Now we proceed to the first LL^ correction to 
the quasielastic scattering of color-triplet quarks. 
As in the case of the free-nucleon target, it must 
be evaluated for Zc = Zg <S^ 1 and Zf, = 1 with 
soft-gluon approximation H87|) for the quark split- 
ting function. Combining together the virtual and 
real corrections we find (the negative valued terms 
originate from VRC) 



— f ^' ' — = /Colog— < S[6;cro(a;o)] j (Pk 2K{p,p - K)(j){b,xo,K) ~ (f>{b,XQ,p)K{K,K - p) 



/ (fKi(fK2(t){b,Xo,K2) 



X K{p+ Ki,p+K2)(j){b,Xo,K.i)- K{Ki,p+Ki+K2)(f>{b,Xo,p) > 

= 5(t>{b,x,p). (104) 



In conjunction with the Born contribution H29|l . 
we obtain the LL^ evolved quasielastic scattering 
spectrum 



d<7Qei{x,p) 
cPbcPp 



ib,x,p), 



(105) 



which is linear factorizable in terms of the collec- 
tive nuclear unintegrated glue as defined through 
the excitation of coherent diffractive dijets. By 
hindsight, we attribute this finding to the Abelian- 
ization of the single-jet problem to the LL^ ap- 
proximation, for the related discussion see Ref. 

m 



D. The first LL^ iteration for nuclear glue: 

octet-octet color dipoles and color-octet 

leading quarks (gluons) 



r 



inclusive spectra of leading color-octet partons off 
nuclear targets. In this case there is no need to in- 
voke the large- iVc approximation. One defines the 
octet-octet (gluon-gluon) collective nuclear glue 
^gg {b,x,p) as the Fourier transform of the nuclear 
S-matrix for octet-octet color dipoles. The Born 
approximation, ^gg{b, xo,p), is given by Eq. H73|l . 
We also recall the convolution property [l7| 



<^gg{b,xo,p) 



$n 



.$J(xo,p) (106) 



Still more insight into the properties of inclu- 
sive production of leading partons comes from the to hold to the Born approximation. 

I 



The first LL^ iteration for the octet-octet nuclear S-matrix reads 

X 

^S,,(6;.,r)./rf.,/rfVl*..(.„p)-*..(.„P + r)|^{s[6;.,,,(a:o)]-S[b;a,,(.o)]} 
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= dZg (fp\-^agiZg,p) -'^ag{Zg,p + r)\^ 

X |S[5; ag{xo, p)]S[b; ag{xo,p + r)]S[b; ag{xa, r)] - S[b; agg{xo,r)]j. 
Proceeding the now familiar route, one would find 



S<^gg{b,X,p) ^2}C0lOg^ 



(107) 



X <^ / (fKsd^Kid'^K2<^g{b, Xo,p~ K3)$g(6, Xq, Ki)$g(6, Xq, K2)K{k3 + Ki,Ks + K2) 
- / d'^K.3d^Kid'^K2^gib,Xo,p- K3)^g{b, Xq, K.3)^g{b, Xq, H2)K {ki, K3 + Ki + K2) > 

= /Colog — < -S^[b;crg,o(a;o)]'5^^^(p) / d^Kid'^K,2(j)gib,xo,Ki)K{K2,Ki + K2) 

r r 

+S'^[b;ag^o{xa)] / d^n 2K{p,p + K)(l)g{b,xo,K) - (l3g{b,xo,p)K{K,p + k.) 
+S [6; cTg, o(a;o)] / d'^Kid^H2 (l)g{b,xa,K.i)(j)g{b,xo,K2)K{p + Ki,p + K2) 

+ 2<j)g{b, Xo, Ki)(j)g{b, Xq, K2)K{p + Ki,p + /ti + K2) - (t>g{b, Xf),p)(j)g{b, Xq, K2)K{ki,Ki + K2) 
-cj)g{b,Xo,p- K2)cj)g{b,X(i,li2)K{Ki,Ki + K2) 



d^K3d^Kld^K2 



,(b,Xo,p- H3)<l>g{b,Xo,Ki)(j)g[b,Xii,H2)K[K3 + /ti,K3 + 1^2) 



/ d^qd^Hid^K24>g{b,Xf),p- K3)4>g{b,XQ,K3)(j)g{b,Xf),K.2)K{Ki,K3 + Ki+ K2) > 



= SSgg^Q (6, X)6'^ (p) + 54>gg {b,X,p) 



(108) 



r 



The LL- evolution of the the triplet-antitriplet 
collective nuclear glue <^{b,x,p),(f){b,x,p) is a 
(quadratic) nonlinear fcx-factorization quadra- 
ture (|100|l . (|102|) in terms of the same quan- 
tity. In striking contrast to that, for octet-octet 
color dipoles 5^gg{b^x.p) is a (cubic) nonlin- 
ear fcx-factorization quadrature of another quan- 
tity - ^g{b,x,p). The nuclear absorption-non- 
suppressed component of pU8ll can be viewed 
as a fusion of three nuclear pomerons described 
by ^g(b,x,p) to a different nuclear pomeron de- 



scribed by ^gg{b,x,p) 

The above distinction between the pomeron fu- 
sion properties of the nonlinear BFKL evolution 
for the triplet-antitriplet and octet-octet dipoles 
was not discussed before. Now we notice an 
important unifying aspect of the LL- evolution 
properties of $(6,a;,p) and ^gg{b,x,p). Namely, 
the same answer for S^gg{b, x,p) could have been 
derived from the already available results for the 
triplet-antitriplet color dipole. Indeed, notice the 
factorization property 



S[6; ag{xo,p)]S[b; <7g{xo,p + r)]S[b; (Jg{xQ, r)] - S[6; agg{xo,r)] = 
= S[b; ag{xQ,r)]i^5[b; ag{xo, p)]S[b- ag{xo, p + r)] - S[b; crg(a;o, r)]|. 

I 



(109) 



The quantity in curly braces in the last line of to a simple substitution 

H1U9|I is identical to the integrand of l|^ subject 

S[b;a{c)\ => b[b;(7g(c)\. 



(110) 
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Let 5^g{b,x,p) be the result of the iteration H98|) 
subject to the substitution H110|l . Then the factor- 
ization property (|109|) amounts to the convolution 
representation 

S<i>gg{b,X,p) = 2(^<i>g{b,X0)(g)d'^g{b,x)j{p)- 

(111) 

The factor 2 in the rhs of Eqs. ((TnS|l . l(TTTl) comes 
from the ratio of Casimirs in the q -^ qg and g -^ 
gg splitting functions. Then Hlll() suggests that 
the evolved ^gg{b, x,p) is a self-convolution IjlUtil) 
of the evolved <&g(b, x,p): 

^gg{b,X,p) = ^gg{b,X,p) + 6^gg{b,X,p) 
= l[^g{b,X0)+6<^g{b,x)] 



i['^g{b,X0)+6<^g{b,x)] ]{P) 



$„ 



$g ]{b,X,p). 



(112) 



Recall that to the leading order of large- TVc pertur- 
bation theory ^g{b,x,p) = ^{b,x,p). Then Eqs. 
(|100() . H102() in conjunction with the convolution 
H112I) provide a unified description of the LL^ evo- 
lution for both the triplet-triplet and octet-octet 
collective nuclear glue. 

A short comment on SSggfi{b) is in order. Be- 
cause S [b;agfi{xo)] = S[b;agg,o{xo)], it can be 
cast in the form 

SSgg^o{b) = -2S[6;crgg,o(2^o)]^olog — 

X / (pK,iCpK2(j)g{b, Xq, Ki)K{k2,Ki + K.2) ■ 

(113) 

Apart from the factor 2, the origin of which has 
been explained following Eq. (|lll|l , it is an exact 
counterpart of Eqs. (PS|l and (|101|l . 

We leave it as an exercise to check that when 
put together, the LL^ approximations for VRC, 
Eq. ITHl, and REC, Eq. (|73), give precisely 



(Pb(Pp 
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{b,x,p). 



(114) 



We again recover the linear fc^-factorization for 
the quasielastic scattering spectrum in terms of a 
judiciously chosen collective nuclear glue. Basi- 
cally, this property is made obvious by Eqs. \iVl\ 
and ITUSJ. 



E. Going beyond the first iteration of the 
LLi evolution? 



Will the exciting finding of linear fc^- 
factorization for quasielastic scattering of partons 
off nuclei hold to all orders of the LL- evolution? 
Will the leading jets offer a long sought linear 
mapping of the collective nuclear glue? 

As far as further iterations of Eqs. H98|l and 
p()7|l are concerned the answer is a negative one - 
such a BK approximation can not be treated as a 
closed form equation. Still here we present heuris- 
tic arguments which partially vindicate the BK 
approximation - our scrutiny of all steps behind 
the derivation of Eqs. (|98|l and (|107|l suggests, 
that a limited number of iterations of these equa- 
tions will be a good working approximation capa- 
ble of describing evolution effects in the limited 
range of energies attainable at RHIC and LHC. 
Our line of reasoning is as follows. 

We consider the scattering of a thin projectile 
(hadron) off a thick nucleus. The thickness of the 
target nucleus is measured by a large parameter 

j/^(b,xo) = ifTo(a:o)r(6), (115) 

which is a number of nucleons in the tube of nu- 
clear matter of cross section ^cro(a;o) at an impact 
parameter b. For instance, realistic models for 
the color dipole cross section suggest ao{xo) ~ 60 
mb |4g and for central collisions with a gold nu- 
cleus one finds UAuib = 0,xo) ~ 6. We treat 
nuclear effects to all orders in this large param- 
eter VAib^xo). The coherent nuclear effects start 
at xq = xa- At this boundary the nucleus can 
be treated as a dilute gas of color-singlet nucleons 
and the resummation of nuclear effects boils down 
to the Glauber-Gribov exponentials for multipar- 
ton S-matrices. 

Recall that inelastic interactions with the free- 
nucleon target are described to the single-gluon 
exchange in the i-channel. The \Aj- evolution is 
treated as an effect of ag Fock states in the projec- 
tile parton a or aag Fock states in the color dipole 
aa. Here emerges sort of a renormalization group 
(RG): upon the integration over gluon variables, 
one would cast the interaction of the three-parton 
Fock state aag in the form of an interaction of 
the aa dipole. The sole effect of integrating out 
the gluon is a renormalization of the color-dipole 
cross section, which is equivalent to a renormal- 
ization of the gluon-density of the target. I.e., at 
each step of the LL^ evolution, the gluon from 
the aag Fock state of the beam is RG reshuffled 
to become a part of the target, while the beam is 
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always described by the lowest Fock state oo. The 
result is the closed-form linear BFKL evolution in 
either the color dipole |33,|33 or the unintegrated 
glue [a, 01 representations reviewed in Sec. V.B. 
Of course, to the single-gluon t-channel exchange 
for inelastic processes, this renormalization pro- 
cedure is a manifestly beam-target symmetric one 
|4lt48]. 

In our discussion of the LL^ evolution for nu- 
clear targets we followed exactly the same RG ar- 
guments: we reabsorbed the effect of the aag Fock 
state of the beam into the renormalization of the 
nuclear color-dipole S-matrix and the correspond- 
ing evolution of the collective glue $00(5,0;, p) of 
the target nucleus. The color-dipole S-matrix and 
its Fourier transform - the collective glue - are well 
defined to all orders in LL^ and are perfectly ex- 
perimentally accessible: the frequently mentioned 
amplitude of excitation of coherent hard diffrac- 
tive dijets ttA -> q{p)q{-p)A equals pi l2ll2^ 
(we cite it for forward excitation) 



M{x,Zq,p) 



X exp(-ipr) [1 - Sqq{b, x, r)] 

-•^■^{Zq,K)^qq{b,X,p- k) 



I / (fK^^{Zq,K) \ 
■ / (fb(j)qq{b,X,p), 



(116) 



where the last approximation holds for jets with 
the transverse momentum much larger than the 
intrinsic momentum of quarks in the pion. Such 
a direct measurement of 4>qq(b,x,p) is possible at 
arbitrarily small values of x = M^^/Wpj^; a use- 
ful review of the current experimental situation is 
found in [43 ■ In principle, the octet-octet collec- 
tive glue (j)gg{b,x,p) is equally measurable in the 
diffractive breakup of glueballs. 

Then, what makes Eqs. ()98|1 and p07|l suspect 
from the RG point of view? What we do is the 
adding of one more radiated gluon to the final 
state and an evaluation of the matching radiative 
correction to the lower order cross section. In the 
operator form, the master formula 14()|l . the iden- 
tification (|41|l of the multiparton S-matrices, the 
whole discussion of the unitarity driven renormal- 
ization H5Q(I . (|98|l and (|l()7|l of elastic S-matrices 
hold at an arbitrary x < xa- The only caveat 
is the Glauber-Gribov factorization for nuclear 



matrix elements of the multiparton S-matrices. 
Specifically, in all our derivations we made an in- 
discriminate use of the property known to all prac- 
titioners in the Glauber-Gribov multiple scatter- 
ing theory ^, ^J , 

n n 

{A\ n S(6, a,)\A) = X{{A\S{b, a,)\A). (117) 
1=1 1=1 

It derives from (i) the nuclear S-matrix being a 
product of the free-nucleon ones. 



SA{b) = \{SN{b-b,), 



(118) 



where bj is the impact parameter of bound nu- 
cleon J, and (ii) the dilute uncorrelated gas ap- 
proximation for the nucleus. 

Each RG reshuffling of the s-channel gluon of 
the aag Fock state from the beam to a parton of 
the target nucleus introduces in the tube of z^a(5) 
nucleons a gluon which interacts with all nucleons 
of the tube and breaks the assumption of the un- 
correlated gas. Specifically, the Glauber-Gribov 
factorization 

S[b, cr5gg(a;o,p,r)] == S[6; cr(a;o,p)] 

X S[b;(j{xQ,p + r)], 

S[^, crggg(a;o,p,'r)] = S[b;ag{xQ,r)S[b;ag{xo,p)] 
X S[b](Tg{xQ,p + r)], 

(119) 

holds only at the boundary xq — xa, i-e., for 
one s-channel gluon, Ug = 1. It will break down 
at smaller x, beyond the first iteration, i.e., for 
rig > 2. While it prevents one from treating the 
BK approximation as a closed form equation for 
the free-nucleon target (for the related recent dis- 
cussion see Ref . J2^ |43 and references therein) , 
arguably, the breaking effects are small if the mul- 
tiplicity of RG gluons is small compared to the 
number I'Aib) nucleons of in the tube of cross sec- 
tion ^ao{x). To be more precise, the number Ug 
of iterations of the BK approximation must be 
bounded from above: 



1 < VAib). 



(120) 



One can quantify Ug from the experimentally 
observed small-x growth of the unintegrated glue 
in the nucleons and/or the DIS structure function. 
In the fully developed BFKL asymptotics 



F2p{x,Q'')^ [ i 



= exp(Aipy), (121) 
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where Y = log - and Ajp is the intercept of the 
QCD pomeron. The expansion 

exp(Aipy) =y^ ( AipF)"- (122) 

rig! 

is an expansion over the multiplicity Ug of s- 
channel gluons. The average multiplicity equals 
(ng) — AjpY and the average rapidity gap be- 
tween gluons equals AY = 1/Aip. Increasing 
the total rapidity Y by AY amounts to increas- 
ing (rig) by unity and, simultaneously, to the in- 
crease of the structure function by the factor e^. 
The experimental data on A]p from HERA ex- 
hibit a steady rise of Ajp from ~ 0.1 at small Q^ 
to - 0.2 at Q2 ^ lOGeV ^ to - 0.3 at Q^ - 100 
GeV ^ 113. For heavy nuclei xa ^ 10~^, which 
corresponds to the nucleus fragmentation rapidity 
span rjA ~ log — ~ 4.5. If Y is the total rapidity 
span, then 



I ^ AjpiY ~ 7u) - I. 



(123) 



At LHC, for jets with p± - 10 GcV the total 
kinematical rapidity span Y ^ 15, which entails 
an estimate Ug ~1 < 2. For mini-jets with p± ^ 3 
GeV Eq. H123|) gives stiU smaller Ug - 1 < 1.5. 
Comparing that to the above estimate for i^Ain 
we conclude that several first iterations of the BK 
approximation would reproduce the gross features 
of nuclear effects in pA collisions in the energy 
range from RHIC to LHC. 

The practical application to radiative stopping 
of ( or the transverse momentum dependent LPM 
effect for ) leading jets in pA collisions at LHC 
would be as follows. Order the s-channel gluons 
gi-.-Qn in energy in the nucleus rest frame: Zi < 
Z2... < Zn- The stopping will be dominated by the 
hard gluon gn , the remaining n — 1 gluons carrying 
small energy would have a negligible impact on 
the energy loss and their principal effect would be 
the LLi evolution of the collective nuclear glue. 
Then the so-evolved collective nuclear glue must 
be used in the formalism of Sec. IV as an input 
for the quantitative evaluation of the stopping of 
leading quark jets caused by radiation of the hard 
gluon Qn. 

Our discussion was limited to a scattering of a 
dilute and small projectile on a dense and large 
target. Still it has a broad range of applicability 
at energies of RHIC and mid-rapidity production 
at LHC. Indeed, in this energy range proton pro- 
ton interactions which are still far away from the 
strong absorption regime of Oei/crtot = 1/2 and a 
single-pomeron exchange for scattering of a dilute 
projectile on a dilute target is a reasonable start- 
ing point |5lL l52| . An extension to multipomeron 



exchanges in proton-proton interactions and inclu- 
sion of pomeron loops remain a major challenge in 
the field (for a review see pql. 



F. Extension to the LL- treatment of 
mid-rapidity production 

While our principal motivation has been the in- 
terplay of real and virtual radiative corrections 
in the LPM effect and quenching of leading high- 
p± jets, the extension of the above LL- consid- 
erations to mid-rapidity high-p_L jets in DIS or 
in proton-nucleus collisions at LHC is straight- 
forward. Here to LL^ the familiar linear k 1 - 
factorization would hold on the beam side - the 
discussion of cancellations of spectator interaction 
effects in the inclusive jet cross section integrated 
over the phase space of spectator partons of the 
beam is found in Refs. [l7l l87l |. We also recall 
that to the Born approximation, the nonlinear 
/j^-factorization simplifies, for LLi real soft gluon 
emission processes, to the linear one [13; IS^- To 
the next order in pQCD, one must consider the 
effect of real radiation in the rapidity span [tja, v] 
between the nucleus fragmentation region and the 
mid-rapidity jet, and include the corresponding 
virtual correction to the Born cross section. 

We recall first the Born spectrum at the 
boundary xq = xa for mid-rapidity gluons with 
transverse momentum p, (pseudo)rapidity rjg = 
log{z/x) and z = p^ jxWaN ^ 1 from pQCD sub- 
processes q ^ qg and g ^> gg |l7| 



{2TTf da a 



dr]g(PpcPb 

a — >ag 

Z / d^/t(/)gg(6,XA,/«)|*(2:,p) - *(Z,P- /«)P. 

(124) 

For the transverse momenta above the infrared 
parameter fig, in H58|) we can use 



K{p,p-k) = 



K 



p2(p-K)2' 



(125) 



identify the unintegrated glue in the incident par- 
ton a, 



dGa{z,p- k) 



2aszPga{z) 



1 



d'^p (p — kY 

make use of H29|) and cast (|124|l in the form 



, (126) 



{2TTfdaA 



drigd^pd^b 



a — >ag 



./.=..W(P-P.-») 
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dGa{z,p^) ^ K?_ ^ daQel{xA,l^) 
d^Pa P^ ^^'^ 

P J 



dGa{z,Pa) ^ dGA^gg{b,XA,K) 

The reabsorption of k? into the definition 
'iiras dGA,gg{b,XA,K) 



(127) 



K (t)gg{b,XA,l^) 



Nr. 



d^ bd^ K 



(128) 



brings p27|l to a lucid beam-target symmetric 
form. The remaining factor 1/p^ is famihar from 
the square of the BFKL gluon radiation vertex 
jo, |3| - an absence of nuclear renormalization of 
this vertex is noteworthy. 

Modulo to the above explained factor K^/p^, 
Eq. H127|l is an exact counterpart of the convo- 



lution (|63|l . Upon this identification the differ- 
ential cross section of quasielastic gluon-nucleon 
scattering, the further analysis of radiative cor- 
rections would not be any different from that pre- 
sented above for leading jets. An important sim- 
plification is that the major effect of real radia- 
tion is a slight shift of the rapidity of the high- 
P_L jet towards the nucleus region. In view of 
the approximate boost invariance, this rapidity 
shift can safely be neglected. Then one is allowed 
to evaluate the real and virtual radiative correc- 
tions in the \Aj- approximation. Next we reiter- 
ate the above point that although the real radia- 
tive corrections are highly nonlinear quadratures 
of the collective nuclear glue 17], upon lumping 
together the real and virtual corrections all non- 
linear effects combine to the nonlinear LL- evolu- 
tion of the collective nuclear glue: 4>gg{b, xa, i^) of 
the Born approximation changes to the nonlinear 
LLi-evolved 4>gg{b, x, k.). 
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FIG. 7: (a) The Kancheli-Mueller diagram for the single-gluon spectrum, (b) its pomeron exchange representa- 
tion for mid-rapidity gluon spectrum in proton-nucleon collisions, (c) its explicit form in terms of the emission 
of real gluons off reggeized gluons qr and free-nucleon gluon density representation for the cut pomeron, (d) the 
pomeron exchange representation for mid-rapidity gluon spectrum in proton-nucleus collisions, (c) the explicit 
form of the diagram (d) in terms of the emission of real gluons with transition between reggeized gluon gR and 
effective nuclear reggeized gluon {gg}R,A and free-nucleon and collective nuclear glue,(f) the representation of 
the diagram (d) in terms of two cut nuclear pomerons TPA,g composed of effective nuclear reggeized gluons gR,A- 



The beam-target symmetric linear k± factor- ization form of Eq. I|127|l shows that sort of a 
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LLi BFKL multiregge factorization |a, Q holds 
for the mid-rapidity single-jet spectrum off nu- 
clei. From the viewpoint of the Kancheli-Mueller 
3^3 optical theorem [Ullll of Fig. [7\ this 
form of the single-jet spectrum suggests an in- 
teresting effective reggeon field theory interpreta- 
tion (whether this interpretation and the following 
conjectures will lead to a viable diagram technique 
|53 | or not needs further scrutiny). In the single 
mid-rapidity jet production on nucleons the free- 
nucleon BFKL cut pomerons IP are exchanged 
between the beam (target) and produced gluon 
(Fig. |2|d), the more detailed structure of this dia- 
gram in terms of the radiation of the gluon from 
reggeized gluons gn is shown in Fig. Et. For the 
nuclear target, Fig. \7\i, on the nucleus side one ex- 
changes the cut nuclear pomeron TPa.qq described 
by dGA.gg{b,x,K)/(Pb(pK which resums all mul- 
tipomeron exchanges enhanced by a large thick- 
ness of the nucleus. If we associate with this col- 
lective glue the effective nuclear reggeized gluon 
{gg}R,A, then this Kancheli-Mueller diagram can 
be cast in the same form, Fig. [7^, as in the free 
nucleon case, Fig. \7jp. Our notation {gg}R,A for 
this effective gluon is a reminder that it can be 
viewed as a composite state of two gluons gR^A as- 
sociated with still another member of the family 
of nuclear pomerons, lPA,g, described by the col- 
lective glue dGA.g{b,x,K)/cfb(pK related by Eq. 
(|128|l to $g(6,x, /«). In due turn, the latter can 
be viewed as sort of a coherent composite state 
of the free-nucleon reggeized gluons. Considering 
the distinct, nonlinear vs. linear, LLi evolution 
properties of the free-nucleon and collective nu- 
clear glue, the possibility of such a resummation 
and representation in terms of a single nuclear 
reggeized gluon {gg}B,.A is far from being an obvi- 
ous one. ( From the color dipole perspective it has 
its origin in the cancellation of spectator interac- 
tions and Abelianization of intranuclear evolution 
of dipoles in the single jet problem.) 

The no-nuclear-renormalization of the BFKL 
gluon emission vertex is a still another nontriv- 
ial property: the gBgng and gR{gg}R,A9 vertices 
are identical. Furthermore, since {gg}R,A is a 
composite state IgR.AdR.A), the nuclear Kancheli- 
Mueller diagram of Fig. [Tji is a triple-pomeron di- 
agram IP —>■ TPA.gTPA,g with three cut pomerons 
( Fig- [Zf)i and the corresponding fully cut triple- 
pomeron vertex is not renormalized by nuclear ef- 
fects and is related to the square of the BFKL 
vertex. If we recall that gR^A is by itself a coher- 
ent composite state of reggeized gluons gR, and 
invoke our experience with unitarity rules for di- 
jet and single-jet production 58], then the above 



non-renormalization property would apparently 
extend to a whole family of fully cut multipomeron 
vertices IP -^ (riIP)(r?T,]P) in the single-gluon pro- 
duction problem. A very different derivation of 
non-renormalization of the gluon emission vertex 
in proton-nucleus interaction, based on the QCD 
shock wave scattering approach, is found in |5q . 
although Balitsky considered only the Born ap- 
proximation and did not discuss the relevant in- 
terplay of the real and virtual pQCD radiative cor- 
rections. 

As we cited above, this multiregge, and linear 
/cj^-factorizable, form of the mid-rapidity gluon 
spectrum has already made an appearance in the 
literature ([13)0)1111 and references therein). A 
useful aspect of our analysis is an elucidation of 
the role of the s-channel unitarity in the evalu- 
ation of the virtual pQCD radiative corrections 
for nuclear targets. Our technique gives a closed 
form nonlinear k± factorization quadratures for 
the spectra without invoking reggeon field theory 
and Abramovsky-Gribov-Kancheli (AGK) unitar- 
ity cutting rules 57J. From the jet phenomenol- 
ogy viewpoint the exchange by a nuclear pomeron 
on the nucleus side of the Kancheli-Mueller dia- 
gram entails a full fledged familiar Cronin effect: 
nuclear attenuation for the transverse momenta 
below the nuclear saturation scale, p^ < Q\, an- 
tishadowing for p^ ^ (2 — 3)Q\ and slow ap- 
proach to the impulse approximation at large p^ , 
controlled by a nuclear higher twist correction oc 
Qa/p^- a full discussion of the hot disputed issue 
of AGK unitarity connections between diffractive 
and inelastic processes goes beyond the scope of 
the present communication (the subtleties of nrul- 
tiple cut pomeron exchanges in pQCD were com- 
mented on in |5q and will be reported elsewhere 
[53). Here we only mention that one version 
of the AGK rules discussed in [B^ suggests the 
substitution (j)gg{b,x,K,) — > ^T{b)f{x,K). This 
amounts to the impulse approximation and misses 
the Cronin effect: numerically the departure from 
the impulse approximation is quite substantial 
(see Ref. [I3 and references therein). While our 
interest is in heavy targets, in his AGK discus- 
sion of single gluon spectra Braun focused on the 
transition from 2 to 4 t-channel gluons relevant to 
the triple-pomeron diagram for the double scat- 
tering in the deuteron ([6l|, [gJ) and references 
therein) . As far as the nuclear departure from the 
impulse approximation is concerned his results for 
double scattering are similar to ours for multiple 
scattering for heavy nuclei. In the reggeon field 
theory formalism used by Braun the LL^ virtual 
pQCD corrections are reabsorbed into the Regge 



trajectory of the gluon. We still find our mani- 
festly unitary approach to virtual corrections an 
instructive one, especially in application to the 
evolution of the color density matrix for collec- 
tive nuclear glue. For instance, when discussing 
the unitarity properties of diffraction, one needs 
to split the qqg contribution to diffractive cross 
section into the pQCD radiative correction to the 
low- to moderate-mass qq diffraction and excita- 
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tion of the high-mass qqg states. Here the cor- 
rect result for the high-mass qqg unitarity cut was 
found to be different from what one would expect 
from the reinterpretation of diffraction as deep in- 
elastic scattering off pomerons treated as a par- 
ticle p3^ . Similar conclusions were reached by 
Braun who correctly noticed how applications of 
AGK rules depend on the form of the so-called 
triple-pomeron vertex. f|6 lll62l|'). 




CC 



CC 




FIG. 8: (a) The Kancheli-Mueller diagram for the cc jet spectrum, (b) its explicit form in terms of the emission 
of a CC pair off reggeized gluons gn in proton-nucleon collisions, (d) the same as (b) for proton-nucleus collisions 
in terms of two cut nuclear pomerons IPa.s composed of effective nuclear reggeized gluons qr^a and the vertex 
of emission of a cc in a transition from reggeized gluon gn to two effective nuclear gluons Qr^a- 
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Now, following our early discussion in |l7l |. we 
comment that the above multiregge form for mid- 
rapidity. Fig. [7|3, is quite likely to be an exclu- 
sive property of single gluon production, and the 
multipomeron version. Fig. Ef, is not an arte- 
fact. "'^ Indeed, a jet can well be formed by a fixed- 
mass two-parton or higher multiplicity multipar- 
ton state. As an illustration we look at the uni- 
versality class of mid-rapidity open charm pair cc 
when the transverse momentum much larger than 
the mass of the diparton. The factorization prop- 
erties of mid-rapidity cc production were stud- 
ied in [12. |l3, 15, 17] . The principal finding is 
a full fledged nonlinear fc_L-factorization for mid- 
rapidity cc production. Consequently, in the cor- 



^ Although a direct connection between the nuclear 
pomerons and next-to-leading order BFKL pomeron is 
not an obvious one, here we cite Fadin et al. who in their 
recent study of the next-to-leading order BFKL stated 
that a rigorous proof of the multiregge factorization in 
elementary scattering can be carried through only for 
one-gluon production ;2li . 



responding 3 — > 3 Kancheli-Mueller diagram the 
exchange between the cc jet and the target nu- 
cleus can not be described by the single nuclear 
pomeron. 

It is sufficient to consider the g —> cc excitation 
and Born cross sections to the large Nc approxi- 
mation. Regarding the color properties, at large 
Nc the cc pair is a composite gluon. Let z and 
z = 1 — z be a partition of the cc jet momen- 
tum between the quark and antiquark, p the total 
transverse momentum of the open charm jet, and 
k the transverse momentum of the quark with re- 
spect to the jet, so that 



Mr - = 

z(l-z) 



p 



(129) 



The differential form of Eq. (23) of Ref. ^\ gives 
the free-nucleon spectrum 



2{2TTfdaN{g* 



cc] 



= 2asPggiz)f{p) 



dzdPpcPk —o-qgy 

X [K{zp + k,k)+ K{(1 - z)p - k, -k)], 

(130) 



30 



Here a square of the vertex of emission of the cc 
pair from a reggeized gluon, [cc)gB.gR in Fig. |S|d, 
is ex Pqg{z)\K{z'p + k,k) + K{{1 - z)p - k, -k)]. 
The nuclear cc jet spectrum equals 



2{2TT)^daA{9* -^ QQ) _ 
<Phdz(Pp<Pk 

= 2asPqg{z) j CpKi(fK25{p - Kl - K2) 
X$g(6, XA, Ki)<i'g{b, XA, K2) 

x[K{zp + k — Ki,k) 
+K{{l-z)p-k^Ki,~k)] (131) 



and corresponds to the Kancheli- Mueller diagram 
of Fig. ^ with two cut nuclear pomerons PA,g- 
Then Eq. H131(l gives an explicit form of the quark 
loop contribution to the fully cut triple-pomeron 
vertex, and a comparison with Eq. 11301) sug- 
gests a simple correspondence between this triple- 
pomeron vertex and the square of the {ccjgngn 
vertex. Eq. (|131|l is the Born approximation, 
the discussion of LL^ iterations in Sec. V.D 
suggests that upon the pQCD radiative correc- 
tions ^g{b,XA,K')^g{b,XA,i^2) is substituted by 
the LLi evolved ^g{b,x,K)^g{b,x,K2), although 
this conjecture must be checked. 

Finally, we briefly comment on a possible struc- 
ture of the related Kancheli-Mueller diagrams 
for digluon jets, a more detailed analysis will 
be reported elsewhere [S^- We recall that the 
mid-rapidity open charm production belongs to 
the universality class of nonlinear fc-factorization 
which is free of coherent distortions of the dipole 
wave functions. Also, at large N^, it is a single- 
channel problem. In contrast to that, intranuclear 
evolution of digluons is a full fledged multichan- 
nel non-Abelian problem, the pattern of nonlinear 
A;-factorization depends on the color state of diglu- 
ons and there are coherent distortions of the color 
dipole wave function which change from one uni- 
versality class to another and depend on the depth 
in a nucleus at which digluons in color represen- 
tations higher than the adjoint representation are 
excited jl^, [iJ, uM ■ The corresponding Kancheli- 
Mueller diagrams for digluons will be similar to 
Fig. IHt complemented by uncut pomeron ex- 
changes on both sides of the unitarity cut, which 
describe the nuclear distortions of the correspond- 
ing partially cut multipomcron vertices. In the 
universality class of digluons in higher color rep- 
resentations these partially cut multipomeron ver- 
tices must be averaged over the depth in a nucleus. 



VI. SUMMARY AND CONCLUSIONS 

We reported a derivation of the transverse mo- 
mentum dependent LPM effect for leading jets 
in the beam fragmentation of proton-proton and 
proton-nucleus collisions. In conjunction with the 
LLi evolution for the collective nuclear uninte- 
grated glue, our results offer a basis for a viable 
pQCD phenomenology of the LPM effect - the 
radiation driven quenching ~ of leading jets in pp 
and pA collisions in the finite energy range from 
RHIC to LHC. 

The first novelty of our work is the derivation 
of the virtual radiative correction to the spectrum 
of leading jets. The quenching of forward jets is 
caused by radiation of hard gluons and, by virtue 
of the unitarity relation, the opening of the radi- 
ation channel is followed by the renormalization 
of the radiationless process. Our derivation of 
this renormalization, i.e., of the virtual radiative 
correction, is based on an explicit solution of the 
s-channel unitarity relation. While for the free- 
nucleon target both the real emission and virtual 
radiative corrections are linear A;j_-factorizable, for 
nuclear targets the quantitative description of the 
LPM effect for forward dijets requires the full- 
fledged nonlinear /cj^-factorization. 

The second novelty of our work is an estab- 
lishment of the LPM effect-deconvoluted, i.e., the 
LLi approximation, spectrum of leading partons 
as a unique linear probe of the collective uninte- 
grated nuclear glue as defined through the excita- 
tion of coherent diffractive dijets off nuclei. This 
must be contrasted to the manifestly nonlinear 
fc^-factorization for all dijet observables and the 
LPM effect for single-jet spectra. Although both 
real-emission and virtual radiative corrections to 
the leading jet spectrum are described by nonlin- 
ear fc_L-factorization, our new finding is that this 
nonlinear fcj_-factorization for leading jets exactly 
matches the equally nonlinear fcj^-factorization for 
the LLi evolution of the nuclear glue. The two 
nonlinear evolving observables prove to be propor- 
tional to each other: the leading jet production 
maps linearly the collective nuclear glue. Our fo- 
cus was on the impact of gluons radiated in the 
rapidity space between the observed leading jet 
and the target nucleus, but all our observations 
on the combined effect of the VRC and REG are 
fully applicable to production of mid-rapidity jets 
too. Here our nonlinear k± factorization results 
suggest interesting reggeon field theory properties 
of the Kancheli-Mueller diagrams for multiparton 
jets. 

For the free-nucleon targets such a linear k±- 
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factorization for, and the BFKL evolution prop- 
erty of, the leading jet spectrum are valid to all 
orders of the LL^ BFKL evolution. While the 
above specified definition of the collective nuclear 
glue in terms of coherent diffraction off nuclei is vi- 
able to all orders of the LL^ evolution, our proof 
of the linear fcj^ -factorization for leading jets off 
nuclei is rigorous only to the first iteration of the 
LLi evolution. Even this result is sufficient for the 

X 

quantitative predictions of the quenching of lead- 
ing jets at RHIC because of the not so large energy 
of RHIC. A similar description of the experimental 
data at the much higher energies of LHC requires 
several - two to three - steps of the IAj- evolution. 
We argue that for such a small number of itera- 
tions, the treatment of heavy nuclei as a dilute un- 
correlated gas of nucleons is still viable, and in the 
finite energy range from RHIC to LHC, one can 
perform this evolution in the Balitsky-Kovchegov 
approximation. After such an evolution, our eval- 
uation of the LPM effect for leading jets will be 
fully applicable to a comparison of pp and pA col- 
lisions at LHC. 

Our analysis of the nonlinear fc_L-factorization 
properties of the BK approximation revealed a re- 
markable feature: the evolution equation can be 
cast in the form of fully linear BFKL evolution 
for the collective nuclear glue complemented by 
a nonlinear term which for hard gluons is a pure 
higher twist correction. To our opinion, this new 
form of the evolution is best suited for several it- 
erations needed to bridge the energy range from 
RHIC to LHC. 

We conclude by noting that one process which 
invites an application of the found linear k±^- 
factorization for leading jets is the coherent 
diffractive deep inelastic scattering scattering off 
nucleons and nuclei. The differential cross section 
of forward diffractive DIS is a sum of differential 
cross section of elastic scattering of multiparton 
Fock states of the photon _38,J34j. As such, diffrac- 
tive DIS is an inherently nonlinear process to be 
described by nonlinear A:^-factorization. One can 
insist, though, on the linear color dipole repre- 
sentation of the diffractive cross section similar 
to that of inclusive DIS [33, |63- Such a rep- 



resentation would entail a certain nonlinear k±^- 
factorization form of the unintegrated glue in the 
pomeron. On the other hand, motivated by the 
above discussion, one may try to define the glue 
in the pomeron in terms of the transverse mo- 
mentum spectrum of leading (anti)quark jets in 
the high-mass diffractive system. The recent find- 
ing is that the two definitions of the unintegrated 
glue in the pomeron would differ one from another 
- one must be cautious with the unitarity cut rein- 
terpretation of the color dipole representation for 
diffraction 1631. 
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Appendix: Still more linear form of the 
nonlinear evolution 



After we partly vindicated the utility of the 
BK approximation for finite-energy evolution and 
heavy nuclei, we report here an interesting alter- 
native form of this equation (the preliminary dis- 
cussion was reported elsewhere |31|), which offers 
a fresh view at the distinction between the linear 
and genuinely nonlinear components of the evo- 
lution equation for the triplet-antitriplet dipoles 
(|102|l . An obvious drawback of this form is a 
reference to the explicit Glauber-Gribov form of 
the nuclear attenuation factor. However, mak- 
ing use of the sum rule / (pK^{b, xo,k) = 1, i.e., 
j cpK(j){b,x,K) = \~S[b]aqqQ(xQ)], ouc cau bring 
it into the form. 



d(l){b,x,p) 
51og(l/x) 



— K.BFKL ' 

+ Q[<i>\{b,x,p) , 



b,x,p) 



(A.l) 



where the nonlinear fcj_-factorization (quadratic) 
functional Q[0] is given by 



J 



Q[4>]{b,x,p) = (fq(fK(j){b,x,q)^ K{p + K,p + q)-K{p,K + p)-K{p,q + p) 

X (p{b,x,K) — (p{b,x,p) K{k,k + q +p) — K{k,k + p) >. 



(A.2) 



The first iteration of (jA.ip with x — xq ~ xa in free of S[b; ao(a;o)] and we view it as a better can- 
thc rhs is exact. The so-defined Q[(/)](b, xojp) is 



didate for finite-energy iterations as discussed in 
Sec. V.E. If one would identify (j){b,x,p) with the 
exchange by collective nuclear pomeron, then Q[(j)] 
describes the fusion of two nuclear pomerons, each 
of which is a nonlinear functional of, and sums a 
multiple exchange by, the free-nucleon pomerons. 
In the case of the azimuthally symmetric 
(f)(b,x,p), i.e., the conformal spin zero in the t- 
channel, the nonlinear fc_L-factorization functional 
Q[(j>] takes a remarkably simple form. In this case 
the azimuthal averaging can be performed explic- 
itly and leads to enormous simplifications. For 
instance, when e^ can be neglected compared to 
all the momenta, we make use of 



P 



(k + p) 



K 



pz 



■0{p' 



(A.3) 
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the related result for a finite e^ in the wave func- 
tion H58() is found in |23,|6J|, see also below. Here 
the subscript indicates the azimuthal angle of the 
momentum the averaging is performed over. 

First we notice that for negligible small e^ = /i^ 



= 2 



K{p + K,p + q) - K{p, K + p) - K{p, q + p) 
p K, + p p q + p 



p^ {k + p)"^ p^ 
K, + p q + p 



1 



(k + p)^ {q + pY p^ 



(A.4) 



Upon the averaging over the azimuthal angles of 
K and q we obtain 



K{p + K,p + q) - K{p, K + p) - K{p, q + p) 



q,K 



-^e[K^-p')e{q^-p^ 



p^ 



(A.5) 



In the second piece of (|A.2p we make first a judi- 
cious shift of the integration variable k in the con- 
vergent integral: K{K,K + q + p) — K{k,k + p) ^ 
K(k — p,k + q) — K{k, k + p). Next, we per- 
form the averaging over the azimuthal angles of 
q and the external variable p. There emerges the 
convergent, and vanishing, integral 



(fK 



1 



1 



0. 



K^ (k + qY 
Upon this cancellation, for negligible small e^ 

^K{K-p,K + q)- K{k,k+p)JJ 



(A.6) 



n K + p 



K p 



K^ 



K? {k + pY {k — pY {n + qY 

,e{K'-p^)e{q'-^^) 



pq 

(A.7) 



Putting it all together the quadratic functional 
Q[(/)] assumes the very simple form. 



Qmb,x,p) = -^ J d\d^Kcl){b,x,q)^{b,x,p)d{q^ -p^)e{K^ -p^) 

-2/Co0(b, x,p)J^J d^qcj){b, X, q)0{K^ - p^)0{q^ - K^) (A.. 

-2}Co(j){b,x,p) -^ d'^q(t){b,x,q), 

Jp^ 1^ Jk^ 



2/Co 

p2 



d q(j){b,x,q) 



valid for the transverse momenta larger than the ^g see equations (jA.15l - lA.19|l below, 
infrared regularization /ig, for the effects of finite 
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It is instructive to reformulate it in terms of 
the collective nuclear glue per bound nucleon, 

2710 s{p^^ 



Hb,x,p) 



One readily obtains 

dJ^A{b,x,p) _ d^GA{b,x,p) 
91ogi dlogp'^dlog- 

= K,BFKL®J^A{b,X,p) 



T{b)TA{b,x,p). (A.9) 



47r/Co 



as{p')T{b){p' 



d?q 



—TA{b,x,q) 



+ TA{b,x,p) \ -^ / -^TA[b,x,q)), 

Jp^ 1^ Jk^ r J 

(A.IO) 

The fusion term is of manifestly negative, i.e., 
shadowing, sign. As such, Eq. IJA.2|) is the bet- 
ter candidate for the definition of of the triple- 
pomeron vertex for nuclear pomerons than Eq. 
(|102|) . Such a representation for the fusion term 
is a new result, it is exact for all momenta, both 
above the saturation scale, p^ > Q\, and un- 
der the saturation scale, p^ < Q\, provided that 

Now we shall argue that at large momenta 
above the the saturation scale the fusion term is of 
pure higher twist. We focus on a; = a;o, where one 
can rely upon the saturation, antishadowing and 
hard momentum properties of (f>{b, xo,p) found in 
\Ul llTl |2J| . Here we need the large-momentum 
approximation 

(t){b,xo,p) w -T(b)/(a;o,p) oc -^— 

(A.ll) 

where the exponent 7^2, see Eq. H14II . Making 
the explicit use of this asymptotic behavior, we 
obtain 

r n 2 

— / (fq(l){b,xo,q) w 
P [Jp^ 

^4>(b,x^,p) I -^ I (fq(p{b,xo,q) 
Jp^ '* Jk^ 

^^^^p'4>'{b,xo,p'). (A.12) 

The fusion term takes the form of the boundary 
condition glue times the higher twist factor. 



Qmb,XQ,p)^ 

j—^cj){b,xo,p)-{p <; 



ib,XQ,p)} 



2Tras{p'^) 

(7 - iras{Q\{b, xo))G{xo,Ql{b, xo)) 

^^^h^.cp{b,x,,p% (A.13) 



pZ 



where the scale for the higher twist is set by the 
saturation scale. In terms of the collective glue 
per bound nucleon, this amounts to 



dTA{b,x,p) 

aiogi 
8/Co^3^g(p^)r(b) ^, 

(7 - lYN.p' 



ICbfkl ^J^Aib,XQ,p) 
T\{b,xa,p). 



(A.14) 



In contrast to the exact integral form (jA.lOp 
of the fusion term, the local in the giuon mo- 
mentum form l)A.14|l is an approximation valid 
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for the specific parameterization 



for al- 
ternative forms of the fusion correction see |6(il 
lerf. The further discussion of HA.10|I and 
(|A.14|I goes beyond the scope of this communi- 
cation. We restrict ourselves to the comment 
that the fusion term does not exhaust the nu- 
clear higher twist corrections - we recall that 
also the boundary condition JF4(b, a;o,p^) con- 
tains a substantial and very similar, but positive 
valued, i.e., antishadowing, linear BFKL-evolving 
nuclear higher twist correction proportional to 
TA{b,xo,p')GA{b,xo,Q'')T{b)/p' MMM 
The shadowing and antishadowing nuclear higher 
twist corrections both have their origin in unitar- 
ity, the impact of the antishadowing correction 
was not given a proper discussion earlier. 

For the sake of completeness, we show the inte- 
gral form of the fusion term in Eq. (jA.lOp subject 
to the infrared regularization (jH^l, i.e. in all wave 
functions. K i n, a ) and JCbfkl we put e^ 
Following |26l l64j we can write 



^^l 



(k + p)2 + ^l^gl K p^ + fig 



2 0(p^/^^M^), 

(A.15) 



where 



^(p^-^/^D = ^ 



X 1 + 



9 9 2 

p - 14 " ** 



(p2 — jjfi — kP'Y + 4p2/^| 



(A.16) 



Eqs. HA.Sp and (jA.5|l take the form 
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K {p + K, p + q) - K {p, K, + p) - K {p, q + p] 



2p2 



(p2 + ^,2)2 



[l-e{p^,.^^?M-e{p^q\^,l) 



q,K 



(A.17) 



K{k — p,k + q) — K{k, k + p) 



2k^ 



p,q (k2+pi2)2 



K^p^,^il)[\-e{^^,q\^il)]. (A.is) 



The infrared-regularization modifies the evolution equation ljA.10|) as foUows: 
dTA{b,x,p) 



91ogi 



= ICBFKL'S>TAib,X,p) 
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TA{b,x,p 



as{,p')T{b)< 



(p2 + ;,2)2 



d?K,d?q 

(k2+ ^2)2^4' 






(/^2,p2^ ^J?^[l - e{K\ q\^i'^)]TA{b, X, q) 



(A.19) 
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